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ABSTRACT 

Using a multiple-lens plane algorithm, we study light propagation in inhomogeneous universes, for 43 
different COiJi?- normalized Cold Dark Matter models, with various values of the density parameter ^Iq, 
cosmological constant Aq, Hubble constant Hq, and rms density fluctuation erg. This is the largest cosmolog- 
ical parameter survey ever done in this field. We performed a total of 3,798 experiments, each experiment 
consisting of propagating a square beam of angular size 21.9" x 21.9" composed of 116,281 light rays from 
the observer up to redshift z — 3. These experiments provide statistics of the magnification, shear, and 
multiple imaging of distant sources. The results of these experiments can be compared with observations, 
and eventually help constraining the possible values of the cosmological parameters. Additionally, they pro- 
vide insight into the gravitational Icnsing process and its complex relationship with the various cosmological 
parameters. 

Our main results are the following: (1) The magnification distribution depends mostly upon Aq and as- 
As (Tg increases, the low-tail of the magnification distribution shifts toward lower magnifications, while the 
high-tail is hardly affected. The magnification distribution also becomes wider as Aq increases. This effect is 
particularly large for models with Aq = 0.8. (2) The magnification probability P,„ is almost independent of 
(Tg, for any combination of flo, Aq, ^^o, indicating that Pm does not depend strongly on the amount of large- 
scale structure. (3) The shear distribution, like the magnification distribution, depends mostly upon Aq and 
(Tg. The shear distribution becomes wider with increasing erg and increasing Aq. The similarities between the 
properties of the magnification and shear distributions suggests that both phenomena are caused by weak 
lensing. (4) About 0.3% of sources have multiple images. The double-image probability P2 increases strongly 
with Aq and is independent of ^o, Hq, and erg. (5) The distribution of image separations depends strongly 
upon Aq, and is independent of erg. Summarizing these results, we find that (1) The properties of gravitational 
lensing, both weak and strong, depend much more strongly upon Aq than any other cosmological parameter, 
(2) magnification and shear are examples of weak lensing caused primarily by the distribution of background 
matter, with negligible contribution form galaxies, while multiple images and rings are examples of strong 
lensing, caused by direct interaction with galaxies, with negligible contribution from the background matter. 
Observations of weak lensing can be used to determine the cosmological constant and the density structure 
of the universe, while observations of strong lensing can be used to determine the cosmological constant and 
the internal structure of galaxies and clusters. Gravitational lensing depends much more weakly upon Qa 
and Hq than as and Aq, making a determination of these parameters from observations more difficult. 
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1. INTRODUCTION 



1.1. Importance of Gravitational Lensing in Cosmology 

The evolution of a homogeneous, isotropic, expanding universe comprised of nonrelativistic matter 
can be described in terms of three parameters: the Hubble constant Hq, the density parameter Qq, and 
the cosmological constant Aq.^ Any combination of these three parameters corresponds to one particular 
cosmological model. The large-scale structure of the universe, galaxies, clusters, superclusters, and voids, 
represents the deviations from this overall homogeneity and isotropy. The most conservative assumption is 
that these structures originate from primordial fluctuations that grow with time as a result of gravitational 
instability. The fluctuations originate from a Gaussian random process, and are characterized entirely by a 
density fluctuation power spectrum. Determining the value of the cosmological parameters, and the correct 
model of large-scale structure formation, is the most important challenge of observational and physical 
cosmology. 

Observing the present, nearby universe can only provide partial information about the value of the 
cosmological parameters, and the power spectrum of density fluctuations (see, e.g. Martel 1995). To unam- 
biguously determine the correct cosmological model, it is crucial to observe the universe at high redshift (i.e. 
at large distances), to study its past structure and evolution, as well as its global geometry. Of all known 
redshift-dependent observational tests, gravitational lensing of distant sources is certainly the most powerful 
and most promising one. 

Because sources must be located at cosmological distance for lensing effects to be important,'' most 
early work on gravitational lensing has focused on QSO's, being the farthest objects in the universe that 
were sufficiently luminous to be observable. However recent developments in instrumentation enable us to 
consider other sources at cosmological distances besides QSO's. The Hubble Space Telescope has deepened 
our field of view tremendously. The Hubble Deep Field (Williams et al. 1996) and Southern Hubble Deep 
Field (Williams et al. 1998) are the deepest sets of exposures even taken. Supcrnovae at rcdshifts z = 0.95 
and z = 1.32 have been observed in these fields (Gilliland, Nugent, & Phillips 1999). These are the farthest 
supernovae ever seen, and as new instruments become available, it will be possible to observe supernovae at 
even larger redshifts, possibly up to z = 20 (Marri & Fcrrara 1998). 

High-z Type la supernovae have the potential to be much more useful than QSO's for gravitational 
lens studies. Unlike QSO's, Type la SNe are nearly standard candles (if we exclude the low-luminosity ones, 
which have limited cosmological use anyway). Although there are some variations in luminosity among 
Type la SNe, these variations are well-understood, and luminosities can be corrected (see the review by 
Filippenko & Riess 1998, and references therein). Having standard candles at cosmological distances has 
major implications for gravitational lensing studies. When a distant source is lensed by the intervening 
large-scale structure of the universe, it is usually magnified or demagnified, without multiple imaging. Only 
a small fraction of magnified sources are lensed sufficiently strongly to have multiple images, as we will show 
in this paper. Without multiple imaging (and without the resolution necessary to resolve the actual shape 
of the image) , a magnified source may be mistaken for a source with a different luminosity, unless the source 
is a standard candle. If it is a standard candle, we can estimate its brightness from its known luminosity 
and estimated distance (based on redshift measurements). If the estimated brightness does not match the 
observed one, then the source is magnified, the magnification factor n is given by the ratio of brightnesses, 
and the combination of the factor /x and the redshift z of the source yields information on the nature of the 
structures responsible for the lensing. 

In recent years, gravitational lenses have been used to estimate or put limits on the values of the 
cosmological parameters. These studies have focused on the cosmological constant Aq (Pukugita, Putamase, 
& Kasai 1990; Turner 1990; Kochanek 1992, 1996a, 1996b; Krauss & White 1992; Maoz & Rix 1993; 
Bloomfield Torres & Waga 1996; Im, Griffiths, & Ratnatunga 1997; Chiba & Yoshii 1997, 1999), the density 
parameter ilo (Yoshida & Omote 1992; Martinez-Gonzalez, Sanz, & Cayon 1997), the (Qq, Aq) plane (Asada 
1997; Park & Gott 1997; Falco, Kochanek, & Mufioz 1998), the deceleration parameter go (Wambsganss et 
al. 1997), the difference Jlo~ Aq (Cooray, Quashnock, & Miller 1999), or the Hubble constant Hq (Watanabe, 

^ Additional parameters must be included if the universe contains additional components such as radiation 
or quintessence. 

^ We are not considering here the possibility of microlensing by nearby massive objects (like MACHO 's). 
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Sasaki, & Tomita 1992; Falco ct al. 1997; Kccton & Kochanck 1997; Kundic ct al. 1997). All these studies 
were focussed on one or two particular cosmological parameters. In any such study, some assumption must 
be made about the value of the cosmological parameters that are not being determined. For instance, studies 
focusing on the cosmological constant all assume a flat universe (flo + Aq = 1), while the studies focusing 
on the density parameter all assume a vanishing cosmological constant (Aq = 0). These assumptions are 
motivated more by theory than observations. The flatness of the universe is a requirement of the standard 
inflationary scenario (Guth 1981; Linde 1982; Albrecht & Steinhardt 1982), while a vanishing cosmological 
constant is the only known way to solve the cosmological constant problem (see, e.g. Weinberg 1989). As 
we shall see, several of these assumptions must be reconsidered in the light of recent observations. 

1.2. The Need for a Full Cosmological Parameter Survey 

During the 1980's, there was a strong theoretical prejudice in favor of the Einstein-de Sitter model 
(ilo = 1, Ao = 0). This model was particularly appealing because it satisfied the flatness requirement of 
the standard inflationary scenario, did not require a cosmological constant, and had fewer free parameters 
than any other model. Numerical simulations of large-scale structure formation (Davis et al. 1985; White et 
al. 1987) showed that an Einstein-de Sitter model in which the bulk of the matter is in form of Cold Dark 
Matter (CDM) could satisfy all the observational constrains known at the time. The theoretical prejudice 
for the Einstein-de Sitter model, supported by the numerical simulations of large-scale structure formation, 
lead to the quasi-universal acceptance of the Standard Cold Dark Matter Model as the model of choice during 
the late 1980's. The only real difficulty that this model faced was the so-called age problem. In the standard 
CDM model (or any Einstein-de Sitter model), Hoto = 2/3. Independent measurements of the Hubble 
constant Hq and the age of the universe showed that this constraint could be satisfied only marginally, as 
one needed to select the smallest possible values for both Hq and to. In particular, the value of Hq had to be 
very near 50 kms~^Mpc~^. Still, this difficulty was not regarded as serious enough to discard the standard 
CDM model. 

Since then, there has been a tremendous increase in the number and quality of cosmological observa- 
tions. These observations fall into three categories: low redshift observations (age of the universe, Hubble 
constant, deuterium abundance, baryon fraction in clusters, abundance of clusters, shape of the cluster mass 
function, large-scale velocity field), intermediate-redshift observations (Type la supernovae, gravitational 
lensing, cluster evolution, Lya clouds), and high-redshift observations (cosmic microwave background). For 
details, we refer the reader to the recent reviews by Steigman, Hata, & Feltcn (1999), Martel & Matzner 
(2000), Ross & Harun-or-Rashid (2000), and Wang et al. (2000). These new observations impose numerous 
constrains on the cosmological parameters and the possible models of large-scale structure formation, and 
invariably argue against the standard CDM model. While the existence of cold dark matter is not ruled out, 
and indeed remains the best approach to explain galaxy and large-scale structure formation, the values of 
the cosmological parameters in the standard CDM model are essentially ruled out, as nearly all observations 
now support a universe with < 1. 

This forces us to consider alternatives to the standard CDM model. A low-density CDM model with 
Oo < 1, known as "Open CDM" or OCDM is in much better agreement with most observations, but does 
not satisfy the flatness requirement of inflation. However, the possibility of "open inflation," or inflation 
without the flatness requirement, has been suggested by several authors (Ratra &: Peebles 1994; Bucher, 
Goldhaber, & Turok 1995; Yamamoto, Sasaki, & Tanaka 1995; Linde 1995; Linde & Mezhlumian 1995), 
and is supported by some recent observations (Filippenko & Riess 1998). But if we stick to the flatness 
requirement of the standard inflationary scenario, then we must include in the model an additional smooth 
component to account for the difference between Qq and 1. Numerous candidates have been proposed (see, 
e.g. Fry 1985; Charlton & Turner 1987), but the most promising ones at present are the cosmological constant 
and the quintessence (see Turner 1999 for a recent review). The corresponding models, ACDM and QCDM, 
have Oo + = 1) where Os is the effective density parameter associated with the smooth component. 

In the standard CDM model, the primordial power spectrum of density fluctuations is assumed to have 
the Harrison-Zel'dovich form P{k) oc fc at large scales. This power spectrum can be modified by introducing 
a "tilt." In this Tilted CDM, or TCDM model, the primordial power spectrum P{k) varies as P{k) oc A;" at 
large scales, where the primordial exponent n can differ from unity. With an exponent n of order 0.8 — 0.9, 
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it is possible to reconcile the constrains imposed by the CMB anisotropics and the present abundance of 
clusters with a density parameter ^Iq < 1, without the addition of a cosmological constant or quintessence. 

To discriminate between these various models, we must determine the values of the cosmological pa- 
rameters. The standard CDM model had no free parameter: the density parameter was unity, and the 
Hubble constant was severely constrained by the age problem. These alternative CDM models have many 
free parameters: the density parameter Oq, the cosmological constant Ao (or the effective density parameter 
associated with quintessence), the Hubble constant Hq, and the tilt n of the primordial power spectrum. 
Gravitational lensing of distant sources provides a powerful tool for studying the structure and evolution of 
the large-scale structure in the universe, and eventually constrain the value of these cosmological parameters. 
However, all the previous studies listed in §1.1 focused on only one or two parameters, while assuming some 
particular values for the other parameters. If the correct cosmological model has several free parameters, it 
is clearly impossible to determine the value of these parameters separately. This would be like looking for 
a fugitive in a city, but limiting the search to one street and one avenue only. A full survey of the cosmo- 
logical parameter space is required in order to determine or limit the values of all cosmological parameters 
simultaneously. In this paper, we present the first study of light propagation in inhomogeneous universe that 
surveys the full 4-parameter phase-space formed by flo, Ao, Hq, and n. 

1.3. Objectives 

There are two distinct objectives to this study. The first one is to determine the properties of lensed 
sources, such as their magnification distributions, shear distributions, image separations, and so on; prop- 
erties which depend upon the value of the cosmological parameters. By directly comparing the predictions 
of the simulations with current and future observations of gravitational lens systems, we hope to eventually 
constrain the values of the cosmological parameters, and the possible scenarios of large-scale structures and 
galaxy formation. 

The second objective is more theoretical in nature. By studying the properties of gravitational lenses, 
and how these properties depends upon the values of the cosmological parameters, individually or in com- 
bination, we can gain insight into the phenomenon of light propagation in inhomogeneous universe and its 
relationship to the underlying cosmological model. This objective is more ambitious than the first one: in- 
stead of merely determining which particular model reproduces observations better, we wish to understand 
the reason which favors this particular model. Achieving this second objective requires that we extend the 
parameter survey to regions of the parameter space that are not particularly favored by observations. This 
forced us to consider a record number of cosmological models in this study, and even with 43 models, our 
covering of the 4-dimensional parameter space is patchy at best. 

The various effects of gravitational lensing can be divided into two broad categories: weak lensing 
and strong lensing.^ Weak lensing is caused by the smooth distribution of matter, with moderate density 
contrast, located between the source and the observer. The magnification or demagnification of sources, or 
the shear circular sources having elliptical images arc examples of weak lensing. Weak lensing provides 
an unbiased information about the matter distribution in the universe, as well as the underlying geometry, 
that can be used to constrain cosmological models (Bacon, Refregier, & Ellis 2000; van Waerbeke et al. 
2000; Munshi & Coles 2000; Jain, Seljak, & White 2000). Strong lensing involves direct interaction between 
the beam and large mass concentrations such as galaxies and cluster of galaxies. Strong lensing causes 
spectacular events such as multiple images, giant arcs, and Einstein rings, providing mostly information 
about the density structure of the lens itself (Schneider, Ehlcrs, & Falco 1992, hereafter SEP, chapter 8), 
though some properties of strong lensing can be used to determine the cosmological parameters as well. For 
instance, the time delay between multiple images of a single source can be used to determine the Hubble 
constant. Most previous studies have focused either on weak lensing or strong lensing. One important goal 
of our work is to study the properties of weak and strong lensing simultaneously. To achieve this, we must 
use an algorithm that can resolve cosmic structures over a very large dynamical range in length, from the 
size of cluster and superclustcrs down to the central cores of galaxies. This lead to the development of a new 
version of the miiltipk^ kms-plane algorithm, based on earlier work by Jaroszyiiski (see §2 below). 



Microlensing would be a third category, which we are not considering in this paper. 
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Because of the nature of these objectives, this projec;t involved an amount of effort that is quite sub- 
stantial compared with similar studies that have been performed and published by various authors in recent 
years. We are studying a record number of cosmological models, 43. We have performed a total of 3,798 
ray-tracing experiments.^ In order to study the properties of weak and strong lensing simultaneously, we 
used beams composed of a very large number of light rays, 116,281. Overall, we have simulated the prop- 
agation in inhomogeneous universes of 441,635,238 light rays. We have generated the images of 3,137,675 
extended sources located at redshift z — 3. For each and every one of these images, we have computed the 
magnification and aspect ratio, and whenever we encountered special kinds of images, such as multiple im- 
ages and rings, we studied their properties, such as image separations, brightness ratios, and hole diameters. 
The calculations started on a Cray J90 supercomputer, and later moved to a more powerful Cray SVl super- 
computer. The A^-body simulations used to generate the large-scale structure for the various cosmological 
models (3 simulations for each of the 43 models, for a total of 129 P^M simulations with one quarter million 
particles each) took about 2000 CPU hours, while the ray-tracing experiments took about 600 CPU hours. 

The remainder of this paper is organized as follows: in §2, we describe the numerical algorithm used 
for the simulations. In §3, we describe the cosmological models included in the study, the ray-tracing 
experiments, and the technique used for analyzing the results. In §4, we review the various elements that 
affect gravitational lensing. Results are presented in §5, including the magnification distributions (§5.1), 
the magnification probability (§5.2), the shear distributions (§5.4), the multiplicity of images (§5.5), the 
distribution of image separations (§5.6), and the properties of Einstein Rings (§5.7). In §6, we discuss the 
various approximations implied by the algorithm and their possible effects on the results. Summary and 
conclusions are presented in §7. 

2. THE MULTIPLE LENS-PLANE ALGORITHM 

Our numerical algorithm was described in detail by Premadi, Martel, & Matzner (1998, hereafter Paper 
I), and convergence tests of the algorithm were presented by Martel, Premadi, & Matzner (2000). In this 
section, we give a brief summary of the method, and describe some minor refinements that have been 
introduced into the algorithm since the publication of Paper I. 

A light ray traveling from a distant source to the observer is affected continuously by the distribution 
of matter it encounters along its trajectory. The multiple lens-plane algorithm (SEF, and references therein) 
consists of approximating this continuous effect by a finite number of instantaneous deflections, caused by 
the matter distribution encountered at various locations along the trajectory of the ray. To implement this 
method, we divide the space between the source and the observer into redshift intervals, and project the 
matter inside each interval onto a plane normal to the line of sight, called a lens plane. Every lens plane 
deflects the light rays that go through it, and the deflection angles can be computed using geometrical 
optics. We can then follow the evolution of a light beam propagating through the universe, by adding 
successively the contributions of each lens plane to the deflection and deformation of the beam. During the 
past decade, this method has been one of the main tools for studying the properties of gravitational lenses 
located at cosmological distances (Blandford & Nayaran 1986; Blandford & Kochanek 1987; Schneider & 
Weiss 1988a, b; Jaroszyiiski et al. 1990; Jaroszyhski 1991, 1992; Babul & Lee 1991; Bartelmann & Schneider 
1991; Wambsganss, Cen, & Ostriker 1996; Bartelmann et al. 1998; Couchman, Barber, & Thomas 1999; 
van Waerbeke, Bernardeau, & Mellier 1999; Hamana, Martel, & Futamase 2000; Jain et al. 2000; see also 
Kochanek & Apostolakis 1988; Paczyhski & Wambsganss 1989; for an interesting alternative, see Fluke, 
Webster, & Mortlock 1999). 

What usually distinguishes a particular version of the multiple lens-plane algorithm from other versions 
is the method used for representing the surface density on the lens planes. The issue of length resolution 
is critical. The magnification of the images of distant sources depends essentially on the amount of matter 
located near the beam, along the line of sight. However, the deformation (or shear) of the images results 
primarily from the tidal influence of distant matter. Therefore, to accurately simulate the effect of both 
magnification and shear, we must reproduce the surface density of the lens planes over the largest possible 
dynamical range in length. Our algorithm achieves a very high length resolution by combining numerical 

® An additional 101 experiments were performed afterward, to answer some specific questions raised by 
an anonymous referee, concerning the angular size of the sources. 
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simulations of largc;-scalc structure formation with a Monte-Carlo method for locating galaxies inside these 
structures. This approach was pioneered by Jaroszyhski (1991, 1992). 



2.1. Large-Scale Structure Formation 



We use a Particle-Particle/Particle-Mcsh (P^M) code (Hockney & Eastwood 1981) to simulate the 
formation and evolution of large-scale structure in the universe. The algorithm produces snapshots of the 
large-scale structure at various redshifts. If we interpret these redshifts in terms of distances from the 
observer, we can treat each snapshot as representing a different region of the universe, and by combining 
them, we can build a chain of cubic boxes representing the large-scale structure over distances of Gigaparsecs, 
from the observer to distant sources. We can then project the matter distribution inside each box onto one 
lens plane. Since different boxes represent different regions of the imiverse, the large-scale structure inside 
neighboring boxes should be uncorrelated. This is clearly a problem if the boxes originate from one single 
simulation, since they would then represent the same large-scale structure at various evolutionary stages. 
To solve this problem, we perform three independent calculations for each cosmological model, by using 
three different sets of initial conditions. We then combine the results of these simulations, such that the 
first simulation provides boxes 1, 4, 7, . . . along the line of sight, the second simulation provides boxes 2, 5, 
8, . . ., and the third simulation provides boxes 3, 6, 9, . . ., thus ensuring that two consecutive boxes never 
come from the same calculation. To reduce correlations even more, we make use of the periodic boundary 
conditions of the simulations, by giving to the matter distribution in each box a random shift. 

For all P'^M simulations, we use 64'' particles and a 128'' grid, inside a cubic box of comoving size 
Lhox — 128 Mpc. The comoving softening length of the algorithm is 300 kpc. The total mass of the system 



is Msys = 3H^noLl^j8nG = 5.821 x 10i'^(f2o/i^)MQ. The mass per particle is m = Msys/643 = 2.220 x 
lO^^(Oo/i^)M0. We performed a total of 129 simulations (3 per model for 43 models). ^° 



The large-scale structure simulations described in §2.1 can be used to compute the effect of distant matter 
on the propagation of the beam. However, at distances less than a few megaparsecs, we cannot ignore the fact 
that matter has collapsed to form galactic-size objects which are much smaller than the resolution of the P^M 
algorithm. We cannot extend the resolution of the P^M code down to galactic scales, because simulating 
the galaxy formation process would require additional physics besides gravity, such as hydrodynamical and 
radiative processes. Instead, we complement the P'^M algorithm with an empirical Monte Carlo method 
for locating galaxies inside the computational volume, based on the underlying distribution of dark matter 
(Jaroszyhski 1991, 1992; Paper I; Martel, Premadi, & Matzner 1998). 

First, we need to determine the number of galaxies present in the computational volume. We assume 
that the present galaxy luminosities follow a Schechter luminosity function, 



where n{L) is the number density of galaxies per unit luminosity. We use the values a = —1.10, = 
0.0156 /i^Mpc"^, and = 1.3 x IO^O/i'^Lq, where h is the Hubble constant in units of lOOkms-^Mpc'^ 
(Efstathiou, Ellis, & Peterson 1988). There is a fourth parameter, the luminosity Lmin of the faintest 
galaxies, which must be introduced to prevent the total number of galaxies from diverging. The value of this 
parameter is not well-known. We assume a value of Lmin = O.OIX* 

Equation (1) allows us to directly compute the present number density no, and luminosity density jo. 



^° These simulations now constitute the core of the Texas P^M Database (Martel & Matzner 2000). 
The universe also contains virialized objects at larger scales, but these are properly simulated by the 
P^M algorithm. 



2.2. 



The Galaxies Distributions 




(1) 




(2) 



(3) 
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where x = L/L^, aJmin = ^min/i*, and T is the incomplete Gamma function. For = 0.01, equations (2) 
and (3) give no = 0.0808 /i^Mpc"^ and jo = 2.13 x lO^/iL© Mpc~^. The total number of galaxies in the 
computational volume is given by 

= no^box = '^*-^boxr(Q: + 1, a;min) • (4) 

For Lbox = 128 Mpc, equation (4) gives iVgai = 28200, 46500, 71500, and 104000, respectively, for the values 
h = 0.55, 0.65, 0.75, and 0.85 considered in this paper (see §3.1 below). 

We use a Monte-Carlo rejection method for determining the location of the galaxies in the computational 
volume at present (z = 0). We divide the computational volume in cubic cells of size (1 Mpc)"', and locate a 
certain number of galaxies in each cell. That number is chosen from a Gaussian distribution with a standard 
deviation proportional to the total matter density in that cell, determined from the distribution of particles 
in the P^M simulation. -'^^ The proportionality constant is chosen in order to reproduce the number of galaxies 
-^gai given by equation (4). The actual location of each galaxy is chosen to be the center of the cell, plus 
a random offset of order of the cell size, also chosen from a Gaussian distribution. With our particular 
choice of cell size, this method naturally creates compact groups of galaxies, with separations comparable 
to the ones found in the Local Group, for instance. Once the position of each galaxy in the computational 
volume at present is known, we reconstruct the trajectories of galaxies, and determine their locations at any 
redshift, by following the trajectory of the P^M particle nearest to each galaxy (a similar method was used 
in Jaroszyhski 1991, 1992). 

Tests have shown that this method produces realistic galaxy distributions (Paper I). In particular, the 
observed 2-point correlation function of galaxies is well reproduced down to separations of order of the cell 
size, 1 Mpc. 

Each galaxy is modeled by a truncated, non-singular isothermal sphere, whose parameters depend upon 
the galaxy luminosity and morphological type. We adopt the galaxy models described by Jaroszynski (1991, 
1992). The projected surface density of each galaxy is given by 

<j{r) = { 4G(r2 + r2)i/2 ' < ''^^^ ' (5) 

, r > Tmax ; 

where r is the projected distance from the center. The parameters rc, Tmax, and v are the core radius, 
maximum radius, and rotation velocity, respectively, and are given by 



. = .o(£)\ (8) 

where the parameters Tco, ^'maxo, vo, and 7 are given in Table 1 (Chiba & Futamase [1999], used a similar 
approach, with different values for the parameters). We use a Monte-Carlo method to generate for each galaxy 
a luminosity L > Lmin, with a probability P{L) proportional to n{L). We determine the morphological 
type of each galaxy by using the observed morphology- density relation (Dressier 1980; Postman & Geller 
1984). Regions of the sky with high concentration of galaxies contain on average more early- type galaxies 
(ellipticals and SO's) and fewer late-type galaxies (spirals) than regions with lower concentration of galaxies. 
By combining this relation with a Montc-Carlo method, we can ascribe a morphological type to each galaxy. 
This can be quite important. Several authors (e.g. Krauss & White 1992; Kochanek 1996a) have found that 
the lensing effect of galaxies is much more important for ellipticals than for spirals, mostly because ellipticals 



The use of a Gaussian distribution is an improvement over the original algorithm presented in Paper I, 
because it allows for the presence of galaxies in low-density regions. 
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have a smalkir eore radius. Consequently, it is not sufEeient to have a reahstie distribution of galaxies: the 
distribution must be realistic within each morphological type, otherwise, some effects, such as strong double 
lensing by pairs of elliptical galaxies, would be underestimated. 



TABLE 1: Galaxy Parameters 



Type 


To {h ^kpc) 


r-maxo {h ^kpc) 


Vq (kms "'") 


7 


Elliptical 


0.1 


30 


390 


0.250 


SO 


0.1 


30 


357 


0.250 


Spiral 


1.0 


30 


190 


0.381 



By combining the distribution of background matter simulated by the P"^M algorithm with the distri- 
bution and surface densities of galaxies, we are effectively describing the surface density of the lens planes 
over 8 orders of magnitude in length, from the size of the largest superclusters and voids, ~ 100 Mpc, down 
to the core radii of the smallest galaxies, ^ 1 pc. The combination of fully nonlinear large-scale structure 
formation, galaxy distributions that reproduce the observed 2-point correlation function, morphological type 
distributions that reproduce the observed morphology-density relation, and galaxy surface density profiles, 
gives to the matter distribution in our algorithm a level of realism that was not present in any of the previous 
studies. 

Figure 1 illustrates the dynamical range that can be achieved with this method. The top left panel shows 

the large scale structure of the imiverse as simulated by the P^M code, for a box of 128 Megaparsecs in size. 
The dots are the particles used by the P-^M algorithm, and are therefore mass-tracers. This panel shows a 
complex network of clusters, filaments, and voids. We enlarge a region of size 10 Mpc and display it on the 
top right panel. At that scale, we see several irregular clusters. The dots still represent the particles used in 
the P-^M code. We enlarge the central cluster, and display it on the middle right panel. In addition to the 
P^M particles, we also plot the galaxies (large dots), whose locations were determined by the Monte-Carlo 
part of the algorithm. We isolate a group of 5 galaxies in this cluster, and plot them at a smaller scale on 
the middle left panel. This group, which is about 3 times smaller than our own "local group," is composed 
of 3 spiral galaxies, one SO galaxy, and one elliptical galaxy. We enlarge the central galaxy, and display it 
on the bottom left panel. The large circle represents the edge of the dark-matter halo (radius r = rmax)- 
The central dot represents the galactic core. We enlarge this core, which has a radius r = Tc, and display 
it in the bottom right panel. There are 7 orders of magnitudes in length between the diameter of this core 
and the size of the largest structures shown on the top left panel, and some galaxies in the simulation have 
a core which is an order of magnitude smaller than the one represented here. 

3. THE COSMOLOGICAL PARAMETER SURVEY 

3.1. The Cosmological Models 

All models considered in this paper are Tilted Cold Dark Matter models (TCDM) , normalized according 
to the results of the COBE DMR experiment. The density fluctuation power spectrum for these models is 
described in great detail in Bunn & White (1997). The power spectrum at the initial redshift Zi is given by 

/ c \ 

P(fc,^0=27r2(^— j 5lC-\zi,Q)k-T^^y,{k) . (9) 

where c is the speed of light, Hq is the Hubble constant, and C{zi, 0) is the linear growth factor between the 
initial redshift Zi and the present, and Tcdm is the transfer function, given by 

TcT,M{q) = ^^^^+3^^'^ [1 + 3.89g + [IQ.lqf + {^AQqf + {Q.llqf] (10) 
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128 Mpc : Superclusters and voids 10 Mpc : Clusters of galaxies 




1.2.8 kpc : Dwarf spiral galaxy 400 pc : Galactic core 




galactic iialo 



Figure 1: Series of zooms illustrating the dynamical range of the algorithm. Top left panel: Distri- 
bution of background matter projected on a lens plane a,t z = 0; Top right panel: Enlargement of a 
dense region, showing several clusters of galaxies; Middle right panel: A particular cluster of galax- 
ies. Small dots represent P'^M, dark matter particles. Large dots represent actual galaxies; Middle 
left panel: A small group of galaxies inside the cluster, composed of 3 spirals, one SO, and one ellip- 
tical. Bottom left panel: A particular spiral galaxy, with a halo radius rmax given by equation (7); 
Bottom right panel: The core of the spiral galaxy, with a radius rc given by equation (6). 
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(Bardeen et al. 1986), with q defined by 



q=[r^y-"\^,h-)-'Ql,, (11) 

a = a'^^'"''^°a~''^'''''^°^\ (12) 
ai = (46.9no/i')°-'™ [1 + (32.mo/i')-"-'"'] , (13) 
a2 = (12.0Jlo/i')°-^'' [1 + (45.0^1o/i')-°-H , (14) 

(Hu & Sugiyama 1996, eqs. [D-28] and [E-12]), where ileo is the density parameter of the baryonic matter, 
and 02.7 is the temperature of the cosmic microwave background in units of 2.7K. The density fluctuation 
5h at horizon scale is given by 



-5o-0-35-0.191nOo-0.17ft„-ft-0.14n^ 



Ao = 0; 



5h={ . ^ . ^^5^ 

i 1.94 X 10-5jl-0-785-0 051nnog-0.95n-0.169n= ^ Aq = 1 - Oq; 

where n = n — 1. 

This power spectrum is characterized by 6 independent parameters: fioj ^^bo, Aq, Hq, Icmb, and n. 
The normalization of the power spectrum is often described in terms of the rms density fluctuation erg at a 
scale of 8/i~^Mpc, defined by 

al = —j^ P{k)W^{M)edk, (16) 

where i = 8/i~^Mpc, and W is the window function, given by 

3 

W{x) = — (sin X — X cos x) . (17) 

The value of erg is a function of the 6 aforementioned parameters. We invert this relation, treating 
CTg as an independent parameter, and the tilt n as a dependent one. We also set Tcmb = 2.7 K and 
^Bo = O.OlSft.^^ for all models, thus reducing the dimensionality of the parameter-space from 6 to 4.-'^^ The 
independent parameters in this parameter space are therefore ilg, Aq, Hq, and as- 

We survey this parameter space by considering 43 different cosmological models. This constitutes the 
largest parameter survey ever done in this field. The values of the parameters are listed in the first 4 columns 
of Table 2, while the values of the dependent parameter n are listed in the fifth column. The values of as 
were chosen by imposing that n remains in the range [0.7, 1.3]. 

3.2. The Ray- Tracing Experiments 

For each model, we performed numerous ray-tracing experiments. The number of experiments for each 
model is listed in the sixth column of Table 2. In each experiment, we compute the propagation of a beam 
consisting of 341^ = 116,281 light rays forming a square lattice on the image plane. The size of the beam 
is 21.9" X 21.9", and the separation between rays is 21.9"/341 = 0.064". This is significantly smaller than 
the typical size of an emitting region. A source of angular diameter 1" will contain 190 rays (more if the 
source is magnified), enough to resolve details such as multiple images. Wc locate the source plane at a 
redshift zs = 3, which is a reasonable choice. The effect of lensing would be more important for sources at 
larger redshifts, but these sources would be more difficult to observe. By choosing zs = 3, we hope to obtain 
results that can be compared with current observations. For a list of the source redshifts used in similar 
calculations, see Table 3 in Martel et al. (2000). The redshift intervals between the lens planes were chosen 
as in Paper I, §3.2.2. 



Recent results support ^bo = 0.0193/i~^ (Buries & Tytler 1998), a value slightly larger than the one 
we assumed. The difference is too small to affect the power spectrum in any significant way, for any of the 
models considered. 
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TABLE 2: Parameters of the Models 





Ao 


Ho 


C8 


n 




0.2 


0.0 


55 


0.3 


1.2187 


99 


0.2 


0.0 


65 


0.3 


1.0966 


92 


0.2 


0.0 


65 


0.5 


1.3188 


104 


0.2 


0.0 


75 


0.3 


0.9993 


47 


0.2 


0.0 


75 


0.4 


1.1228 


199 


0.2 


0.0 


75 


0.5 


1.2190 


70 


0.2 


0.0 


75 


0.6 


1.2979 


51 


0.2 


0.0 


75 


0.7 


1.3648 


199 


0.2 


0.0 


85 


0.3 


0.9191 


101 


0.2 


0.8 


55 


0.8 


1.2057 


99 


0.2 


0.8 


65 


0.6 


0.9326 


70 


0.2 


0.8 


65 


0.7 


1.0062 


51 


0.2 


0.8 


65 


0.8 


1.0702 


108 


0.2 


0.8 


65 


0.9 


1.1269 


50 


0.2 


0.8 


65 


1.0 


1.1568 


67 


0.2 


0.8 


75 


0.6 


0.8273 


51 


0.2 


0.8 


75 


0.8 


0.9629 


51 


0.2 


0.8 


85 


0.8 


0.8749 


45 


0.5 


0.0 


65 


0.8 


0.9457 


48 


0.5 


0.0 


65 


1.0 


1.0439 


77 


0.5 


0.0 


75 


0.8 


0.8686 


51 


0.5 


0.0 


75 


1.0 


0.9656 


50 


0.5 


0.5 


65 


0.8 


A TO AO 

0.7808 


1 A CT 

105 


0.5 


0.5 


65 


1.0 


0.8807 


62 


0.5 


0.5 


75 


0.8 


0.7049 


51 


0.5 


0.5 


75 


1.0 


0.8024 


51 


0.7 


0.0 


65 


0.9 


0.8461 


96 


0.7 


0.0 


65 


1.1 


0.9346 


60 


0.7 


0.0 


75 


0.9 


0.7773 


50 


0.7 


0.0 


75 


1.1 


0.8648 


50 


0.7 


0.3 


65 


0.9 


0.7720 


64 


0.7 


0.3 


65 


1.1 


0.8601 


77 


0.7 


0.3 


75 


0.9 


0.7042 


51 


0.7 


0.3 


75 


1.1 


0.7912 


46 


1.0 


0.0 


55 


1.0 


0.8465 


145 


1.0 


0.0 


65 


0.9 


0.7234 


184 


1.0 


0.0 


65 


1.0 


0.7698 


156 


1.0 


0.0 


65 


1.1 


0.8120 


175 


1.0 


0.0 


65 


1.2 


0.8506 


168 


1.0 


0.0 


65 


1.3 


0.8861 


176 


1.0 


0.0 


75 


1.0 


0.7094 


142 


1.0 


0.0 


75 


1.2 


0.7893 


58 


1.0 


0.0 


85 


1.0 


0.6605 


51 
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3.3. Analyzing the Distribution of Light Rays 

To analyze the results of the experiments, we lay down on the source plane a square grid composed of 
31 X 31 = 961 square cells. The location and size of the grid are adjusted such that, in the absence of lensing, 
the edges of the grid would correspond to the edges of the beam. Figure 2 shows the beam configuration and 
the grid on the sourc:e plane, for a typical experiment. In the absence of lensing, each cell would contain 121 
rays. A source located in a cell containing more than 121 rays would be magnified, whereas a source located 
in a cell containing less than 121 rays would be demagnified. By counting the number of rays in each cell, 
we can then compute a magnification map. The magnification /ie,i in cell i for experiment e is then given by 



Ne,i 
(AT)' 



(18) 




Figure 2: Configuration of the beam on the source plane, for a typical experiment. Each dot 
represents one light ray. The square grid that we lay down on the source plane is also indicated. In 
the absence of lensing, the beam would coincide with the grid, and each cell would contain 11x11 
light rays. 



where Ne,i is the number of rays in cell i for experiment e, and {N) is the average number of rays per cell. 
Looking at Figure 2, we see many cells located along the edges of the grid that are empty or partially empty. 
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This is of course a immcrical artifact caused by the finite size of the beam. Thc«c; c;clls should be cxchided 
from the analysis, otherwise they would artificially bias the magnification distributions toward low values. 
Our original approach was to exclude all cells located along the edges of the grid, using only the inner 29 x 29 
cells instead of the full 31 x 31. However, this proved insufficient in cases of large beam deformation, such 
as the case shown in Figure 2. For this reason, we exclude all cells that are along the edges of the grid, plus 
any cell that is adjacent to an empty cell, since that cell might be half-empty, or even completely empty. In 
Figure 3, we show an enlargement of the right edge of Figure 2. The cells located on the left of the thick 
line are included in the analysis, but the ones on the right are excluded: 6 of them are along the edge of the 
grid, and 3 others are adjacent to empty cells. 



Figure 3: Enlargement of the right edge of Figure 2. The thick line separates cells that arc included 
in the analysis (on the left) and cells that are excluded (on the right), either because they are along 
the edge of the grid or because they are adjacent to an empty cell. 

The number of rays exceeds the number of cells by a factor of 121. There are two reasons for using 
such a large number of rays. First, if the number of rays were comparable to the number of cells, the 
magnification map would be very sensitive to the actual location of the grid. With an average of 121 rays 
per cell, the sensitivity to grid location is significantly reduced. Second, with a large number of rays per cell, 
we can study the properties of the images. Each cell on the source plane constitutes a potential location for 
a source. Assuming that a particular cell contains a source, we can identify the light rays contained in that 
cell and trace back these rays on the image plane. This enables \is to study the properties of images, such as 
multiplicity, angular separations, brightness ratios, and shapes (simple images, arcs, rings, . . .). An obvious 
problem is that the cells are squares, while sources are expected to be circular. We solve this problem by 
superposing over each square cell a circular cell whose diameter is equal to the diagonal of the square cell, as 
shown in Figure 4. These circular cells have an angular diameter of 1", and in the absence of lensing, each 
one would contain {N) = 121(7r/2) = 190 rays. 

4. THE ELEMENTS OF GRAVITATIONAL LENSING 

Before we discuss the results of the experiments, we first review the various elements that enter into 
gravitational lensing, and the relationship between these elements and the cosmological parameters. This 
will facilitate the interpretation of the results presented in the following section. 
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Figure 4: Dashed lines: portion of the square grid laid down on the source plane. Solid circles: 
Circular cells superposed on the square cells. 



4.1. The Cosmological Distances 

The angular displacement of rays caused by lensing depends critically on the angular diameter distances 
between the source and the observer, Ds, the source and the lens, -DlS) and the lens and the observer, 
D]^. In our simulations, there is a continuous distribution of matter between the source and the observer 
(approximated by a finite number of lens planes). In this context, the term "lens" refers to the fraction of that 
matter distribution that is primarily responsible for lensing. The amount of lensing depends on the product 
-Dl^ls, which is maximum for Z?l — -Dls- Hence, lensing is caused primarily by the matter located about 
half-way between the source and the observer. As we discussed in Paper I, this distance effect is competing 
with two others effects: First, the large-scale structure, which is responsible for lensing, grows with time, 
favoring lenses that are closer to the observer. Second, for a fixed number of galaxies, the number density of 
galaxies decreases with time as the universe expands, and therefore the beam is more likely to hit a galaxy 
at locations that are closer to the source. We found in Paper I that the distance effect dominates over the 
other effects, and therefore the matter located half-way between the source and the observer is responsible 
for most of the lensing. The only exception is the Einstein-de Sitter model, in which the large-scale structure 
keeps growing all the way to the present. The matter responsible for most of the lensing tends to be located 
somewhat closer to the observer than the half-way point. This effect is negligible for other models because 
the large-scale structure does not grow all the way to the present, but instead freezes-out at some redshift 
of order Zh ~ 1/flo- 

If the distance effect dominates, then the source redshift determines the "lens" redshift, and therefore 
only one distance, say Ds, is independent. This distance is given by 

Ds = ^f{no,Xo,0,Zs), (19) 



where the function / is given by: 



2 



(2 - Oo + iloZj){l + ^oZiY^^ - (2 - Oo + fio^i)(l + ^ozj)^ 



/2 



f{no, 1 - no, zi, zj) = I dz Oo(l + zf + (1 - %)] (21) 



1 r- 
1 + L, 



(Fukugita et al. 1992). Figure 5 shows the angular diameter distances for the cosmological models considered 
in this paper. The distances depend on Q,q, Aq, and Hq, but not a^. The dependence upon Aq can be quite 
strong, which is why gravitational lenses can impose very stringent limits on the value of the cosmological 
constant (e.g. Fukugita et al. 1990; Kochanek 1996a). 
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12 3 4 5 

z 

Figure 5: Angular diameter distance Ds in units of the Hubble radius c/Hq, versus redshift z, for 
the various models considered in this paper. Solid curves: matter-dominated models (Aq = 0); 
dotted curves: flat, cosmological constant models (Oq + Ao = 1). Curves are labeled with the values 
of Oo- 



4.2. The Mean Background Density 

The importance of gravitational lensing depends on the mean density of matter between the source and 
the observer, which is proportional to the present mean density of the universe, 

^^^^ 

The mean density depends on Hq and ^^o, but not Aq and erg. 



4.3. The Large-Scale Structure 

If the universe was perfectly homogeneous, there would be no gravitational lensing. It is the presence of 

density inhomogeneities, resulting from the growth of the large-scale structure, that is responsible for lensing. 
The amount of structure in the present universe is characterized by the parameter cxg. However, what is 
relevant is not the present large-scale structure, but the large-scale structure at redshift corresponding to 
an angular diameter distance Dl — Ds/2 where most of the matter responsible for lensing is located. Using 
linear perturbation theory, we can estimate the rms density fluctuation as,h at that redshift, 

where £{zh,0) is the linear growth factor between redshifts zi, and 0, which depends on Qq and Aq. Thus, 
cs.L depends on as, ^o, and Aq, but not Hq. For the Einstcin-de Sitter model, £{z,0) = 1 + z, hence erg 
decreases monotonically with increasing zl- For Oq < 1 models, the perturbation "freezes-out" at redshift 
Zfr ~ l/r^o, and grows very slowly between z^ and the present. Hence, at fixed ug, as,L increases with 
decreasing Qq, and if z^ '> z^, (ts,l should be comparable to ag. 



5. RESULTS 

5.1. The Magnification Distributions 

We compute the magnification distributions by combining all experiments within each cosmological 
model, and binning the cells on the magnification maps according to the value /ie,i of the magnification in 
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each cell. We use magnification bins of equal width A/x = 0.02. The probability P{iJ.)Afi that a source has 
a magnification between /U and /x + A/x is given by 



(24) 



where ne(^i, ^2) is the number of cells in experiment e with magnification between and /i2, and the sums 
are over the experiments. Notice that rig (0, 00) is the number of cells for experiment e that are included in 
the analysis, which can vary among experiments (see discussion of Fig. 3). Equation (24) implies 



where the "integral" is actually a sum over bins. By substituting equation (18) into equation (24), we can 
easily show that 



Equations (25) and (26) express the conservation of probability and the conservation of flux, respectively. 

To estimate the accuracy of the distributions, we used a convergence criterion. The magnification and 
shear distributions were always updated each time new experiments were performed in a given model. We 

1 /2 

also computed, in each magnification bin, the uncertainty on the mean AP = <T/iVexp, where cr is the 
standard deviation in that bin, and A'oxp is the number of experiments included in the average. Eventually, 
we reached a point where adding new experiments did not change the distribution, in which case convergence 
has been achieved. For models with large Oq) the convergence is reached much more slowly, and for a few 
models the especially largo number of experiments done was still not sufficient to reach a tight convergence. 
We refer to models for which convergence has or has not been achieved as having "good statistics" or 
"poor statistics," respectively, and we shall be careful when drawing conclusions for models which have poor 
statistics. In Figures 6 9, we compare the magnification distributions for various models. Error bars indicate 
the uncertainty AP on the mean (for clarity, we display error bars for only 1/4 to 1/3 of the bins). The 
important thing to point out is that the error bars are shorter than the separations between curves, except 
in regions where the curves are nearly identical. The differences between the various curves in Figures 6—9 
are therefore real, and not a consequence of insufficient accuracy (except when stated otherwise). 



Figure 6 shows the magnification distributions for three different models with various combinations of 
fio, Ao, and Hq. The various curves represent various values of o-g. For these models, we actually have 
results for 5 different values of ag, but for clarity we only plot the distribution for the smallest, median, and 
largest values of erg. In the absence of lensing, these distributions would be (5- functions centered at /z = 1. 
These distributions reveal that most sources are slightly demagnified, few sources are strongly magnified, 
and the effect of lensing increases with increasing erg. For the open and A models, we see a clear trend: as erg 
increases, the peak of the distribution decreases, the low edge of the distribution moves to even lower values 
(more demagnification) , but the right edge is hardly affected. We see a similar trend for the Einstein-de Sitter 
model, but more experiments are needed to improve the statistics. The explanation resides in the fact that 
the magnification is caused primarily by the matter located near the beam, whereas the matter located far 
from the beam is primarily responsible for the shear. Each lens plane contains a certain distribution of 
matter, with overdense and underdense regions. If the beam propagates through an underdense region, it 
will diverge, resulting in demagnification, and if it propagates through an overdense region, it will converge, 
resulting in magnification. In a cosmological model like CDM, structure formation proceeds hierarchically. 
Small structures form first, then merge to form bigger structures, which merge into even bigger structures, 
and so on. As time goes on, clusters become more massive and less abundant, and the voids between 
them become larger. A larger ag implies that this hierarchical merging process is more advanced, and this 
affects the matter distribution in two ways: first, the underdense regions become more underdense and 
the overdense regions become more overdense, and second, the fraction of the plane covered by underdense 




(25) 




(26) 



5.1.1. The as Dependence 
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regions (the "filling factor") increases while the fraction covered by overdensc regions decreases. In the case 
of demagnification, these two effects act in the same direction: As as increases, the beam is more likely to 
propagate through an underdense region (because of the larger filling factor), and if it does, it will result 
in stronger demagnification, because these regions are more underdense. In the case of magnification, these 
effects act in opposite directions: as as increases, the beam is less likely to propagate through an overdense 
region, but if it does, the magnification will be stronger. As we see, these two effects are almost perfectly 
cancelling each other, and the distributions at values of /i > 1 are essentially independent of as. at least for 
the open and A models. For the Einstein-de Sitter model, the statistics are poor, even though we performed 
significantly more experiments for this model than the other ones. 
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Figure 6: Magnification distributions for various combinations of Qq, Aq, and Hq, showing the effect 
of varying as- Error bars indicate the 1 — a uncertainty AP on the mean, for a few representative 
bins. See text for details. 

5.1.2. The Hq Dependence 

Figure 7 shows the magnification distributions for three different models with various combinations of 
f2o, Ao, and erg. The various curves represent various values of Ho- For the Einstein-de Sitter model (top 
panel) and A- model (bottom panel), the distributions become narrower as Hq increases (the distributions 
of the top panel are quite noisy, but if one focuses on the ordering of the curves along the left edge of the 
distribution, the trend becomes clear). This trend was expected, since the cosmological distances are shorter 
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in models with largo Hq, resulting in a weakening of the effect of gravitational lensing. Notice that this 
effect is partially compensated by the fact that at fixed Qq, the mean background density increases like Hq. 
Prom the point of view of the algorithm, a model with larger Hq has a smaller number of denser lens planes 
between the source and the observer. For the open model (middle panel), these two effect almost perfectly 
cancel each other, and the distributions are essentially independent of Hq. 




Figure 7: Magnification distributions for various combinations of Clo, Ao, and as, showing the effect 
of varying Hq- Error bars indicate the 1 — a uncertainty AP on the mean, for a few representative 
bins. See text for details. 

5.1.3. The ilo Dependence 

Figure 8 shows the magnification distributions for two different models with various combinations of 
Ao, Hq and as. The top panel shows matter-dominated models (Aq = 0), while the bottom panel shows flat 
models (fio + Aq = !)■ The various curves represent various values of Qq- Of the various dependenc;es, the 
fio dependence is the most difficult one to interpret. The reason is that fl^ is the only parameter which all 
elements of lensing, distance, mean background density, and structures, depend on. We are therefore dealing 
with three concurrent effects. As flo increases, the mean background density increases, favoring stronger 
lensing effects. However, the cosmological distances decrease, favoring weaker lensing effects, and at fixed 
as, the large-scale structure at high-redshift is less developed because freeze-out occurs later, also favoring 
weaker lensing effects. For matter-dominated models (top panel), the importance of lensing increases with 
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Oo: resulting in a shift of the distribution toward lower vahies. The dominant effect in this regime is therefore 
the mean background density. The bottom panel shows the magnification distributions for flat models, with 
Ao increasing as SIq decreases. As Figure 5 shows, at redshifts z <3 the angular diameter distances increase 
with Ao for all models, and this reinforces the dependence upon Qq. The distances get even larger with 
smaller fio, and this helps overcoming the mean background density effect. The bottom panel in Figure 8 
shows that, in the limit of large Aq, the effect of having large distances dominates over competing effects. 
However, for values Qq > 0.5 (or Ao < 0.5), there is no clear dependence of the magnification distribution 
on Qq, and we would need much better statistics to determine whether there is an effect. In any case, the 
effect would probably be quite small. 




.7 .8 .9 1 1.1 1.2 1.3 1.4 1.5 1.6 




Figure 8: Magnification distributions for various combinations of Aq, Hq, and as, showing the effect 
of varying Hq. Error bars indicate the 1 — a uncertainty AP on the mean, for a few representative 
bins. See text for details. 

5.1.4. The Ao Dependence 

Figure 9 shows the magnification distributions for two different models with various combinations of 
flo, Hq, and as- The various curves represent various values of Aq. The presence of a cosmological constant 
increases the effect of magnification by increasing the cosmological distances. For the low-density model 
(rio = 0.2, top panel), the cosmological constant results in a widening of the distribution, and a shift toward 
lower magnifications. The effect is much less important for the fio = 0.5 model (bottom panel). This 
indicates that the effect of the cosmological constant is small for values Aq < 0.5, and becomes important 
only for larger values. As Figure 5 shows, the effect of Ao on the angular diameter distances is relatively 
small for models with > 0.5 (or Aq < 0.5), but the difference between the = 0.2, Aq = 0, and the 
fio = 0.2, Ao = 0.8 models is very large. 

5.1.5. Large-Scale Structure versus Galaxies 

In our simulations, the density inhomogeneities responsible for lensing consists of the large-scale structure 
of the universe and the galaxies that are embedded inside that structure. Including large-scale structure and 
galaxies in the algorithm enables us to study both the effect of weak and strong lensing. While most of 
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the weak Iciiising results from the presenee of the large-scale structure, the rare, high-magnification events 
usually result from the presence of a massive galaxy along the line of sight. One may wonder if strong lensing 
is affected by the cosmological parameters, ag in particular. This parameter measures the amplitude of the 
density fluctuations in the large-scale structure, but since the algorithm chooses galaxy locations according 
to the background density, galaxies are more clustered in models with a larger erg- This does not change the 
probability that a given source will be lensed by a galaxy, but it could increase the probability of a source 
being lensed by two or several galaxies clustered together, resulting in a very large magnification. Lensing 
by multiple galaxies must be invoked to explain the properties of some of the images (see Fig. 13 below). 
However, the high-tail of the magnification distributions shown in Figure 6 appears to be independent of the 
value of (Tg. In particular, we do not find the dramatic increase in the values of /U with higher erg that one 
would expect if sources are lensed by compact groups of galaxies. 




Figure 9: Magnification distributions for various combinations of CIq, Hq, and ag, showing the effect 
of varying Aq . Error bars indicate the 1 — a uncertainty AP on the mean, for a few representative 
bins. See text for details. 

We investigate this question by estimating the likelihood that a source will be lensed by two galaxies on 
the same lens plane. Our algorithm models galaxies as truncated isothermal spheres with a mass-dependent 
radius rmax given by equation (7). We computed the geometric cross section Trr^^^ fo^' galaxies, and 
added them up. to find out what fraction /gai of the lens plane is covered by galaxies. For lens planes located 
near z = 0, we obtain 0.00326 < /gai < 0.00503 (the value depends on Hq] see comment following eq. [4]). 
Hence, only a small fraction of the plane is covered by galaxies. However, for sources at = 3, the lensing 
is caused primarily by planes located near z = 1. Since at z = 1 the planes are smaller but the galaxies 
have the same size,^^ we gain a factor (1 -|- z)^ = 4, and the fraction covered by galaxies is then in the 
range 0.013 < /gai < 0.020, still quite small. Galaxy clustering will create regions where the covering is 
larger. If we assume that clustering could increase the space density of galaxies in some regions by a factor 
of 1000 (quite optimistic at z = 1), the surface density would go up by 100, and there could be some galaxies 
overlapping. 

Our current algorithm neglects galaxy evolution, and assumes that the galaxy parameters given by 
equations. (6)-(8) are independent of redshift. 
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Figure 10: Top panel: Galaxy distribution at z = 0.964 for flo = 0.2, \q = 0.8, Hq = 65 kms-^Mpc-i, 
as = 0.8 model. Each galaxy is represented by a dot. Middle panel: enlargement of a galaxy-rich 
region. Bottom panel: enlargement of a rich cluster of galaxies. Galax;ies in the middle and bottom 
panels are represented by circles of radius rmax, as defined by equation (7). 

Figure 10 shows the galaxy distribution on a lens plane located at ^; = 0.964, for the model Oq = 0.2, 
Ao = 0.8, Hq = 65 kms~"'^Mpc~"'^, erg = 0.8. The top panel shows the galaxy distribution over the entire 
lens plane, with galaxies represented as dots. Because of freeze-out, the galaxies at that redshift are almost 
as clustered as at redshift z = 0. We zoom-in on one of the densest, most galaxy-rich region, and display 
that region in the middle panel of Figure 10, where galaxies are now represented as circles with radii equal 
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to Tmax- Becausc the galaxies follow a Schcchtcr distribiition. most galaxies have small radii and few have 
large radii. We zoom-in on the richest galaxy cluster in this region, and display it in the bottom panel of 
Figure 10. There are some galaxies overlapping, but not many. The beam, displayed next to the bottom 
panel of Figure 10, is larger than the galaxies at that redshift, so it might hit several galaxies, but the light 
coming for one individual cell in the beam is unlikely to hit more than one galaxy, in spite of the fact that 
galaxies are strongly clustered in this region. Of course, we have considered only one plane. A light ray can 
hit a galaxy on one plane, and then another galaxy on another plane. But since the structures on neighboring 
planes are uncorrelated, the probability of occurrence of such double hit is unaffected by clustering. The high 
tail of the magnification distributions depends on the distribution of galaxy properties (masses, truncation 
radii, core radii), but not on their actual locations or level of clustering. 



In a large, representative region of the sky, P„i represents the fraction of sources that are magnified. In 
Figure 11, we plot Pm vs. as, for all models. We do not find any particular trend. Instead, Pm is essentially 
independent of erg. This could have been anticipated from Figure 6, which showed that for most models, 
P(/i) is independent of (Tg in the integration range of equation (27). However, this result is in conflict with 
the argument we presented in §5.1.1. We argued that as as increases, the filling factor of overdense regions 
decreases, so a particular source is less likely to be magnified, but if it is, the magnification will be larger. 
However, Pm measures the fraction of sources that are magnified, and does not depend upon how large the 
magnifications are. Hence, according to this argument, Pm should simply be a measure of the filling factor 
of overdense regions, and therefore should decrease with increasing as at fixed Oq, Aq, and Hq. We do not 
see this trend in Figure 11. We even see the opposite trend (Pm increasing with as) in a few cases, such as 
flo = 1, Ao = 0, Hq = 65kms~^Mpc~^ (open triangles in bottom panel). This shows that the interpretation 
of the magnification distributions in terms of filling factors of overdense and underdense regions might be 
sufficient to explain the properties of P(/u) for /x < 1, as we did in §5.1.1, but is too simplistic in the regime 



The basic flaw in this argument is the assumption that Pm depends only on whether or not sources are 
magnified (by having the beam going through a high-density region), and not by how much. If there was 
only one lens plane between the source and the observer, this argument would be correct: Pm would be equal 
to the filling factor of the overdense regions on that plane. However, there are many planes, and even though 
there is an optimal redshift zl where most of the matter responsible for lensing is located, there should be 
several planes at redshift z ^ zl that contribute to lensing significantly. Suppose that the light coming from 
a source hits a high-density region on plane i and then a low-density region on plane j. Plane i will magnify 
the source and plane j will demagnify it. Whether the demagnification caused by plane j is sufficient to 
bring the value of fx below unity will depend on how large the magnification by plane i was. For small values 
of as, the filling factor of overdense regions on plane i is large, and many sources are magnified, but they are 
only slightly magnified (/i ^ 1), and the demagnification caused by plane j will bring many of these sources 
down to fj, < 1. With larger as, fewer sources are magnified, but they are magnified to larger values of fj,, 
and more of them can "survive" the demagnification caused by plane j. This can also be regarded as a 
consequence of the absence of correlation between planes. Iiic;reasing as increases the density contrast on 
each plane, but by stacking uncorrelated planes along the line of sight, the resulting mean density contrast 
is much smoother than the density contrast on any individual plane, and the effect of increasing as is mostly 
averaged-out. 



5.2. The Magnification Probability 



The magnification probability is defined as 




(27) 



/Lt > 1. 
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Figure 11: Magnification probability Pm vs. as- The values of fio and Ao are indicated in each 
panel. The various symbols correspond to various values of Hq. 

The only significant trend we found is that increases with Hq, as shown in Figure 12. In spite of 
the large scatter, the trend is significant. For Qq = 0.2, P„ is between 0.2 and 0.3, while for Oq = 1) Pm 
is larger than 0.29 for all models, with a mean value of order 0.36. The dependence of the distances upon 
Qq is probably responsible for this effect. As Hq goes up, the distances become shorter, and correspondingly 
there are fewer lens planes near 2; ~ 1, reducing the "averaging-out" of the density contrast discussed above. 

Figure 12 shows that for fl < 0.5, Pm does not seem to depend upon Aq. For CIq = 0.7, however, P^ is 
significantly larger for Aq = models than Aq > models. This is very strange, since the Aq = models and 
flat models resemble each other more closely as Hq increases. This effect is probably spurious, a consequence 
of poor statistics. As JIq increases, more experiments must be performed in order to obtain good statistics. 
This was our motivation for performing a large number of experiments for the Einstein-de Sitter model. 
However, we did not performed a particularly large number of experiments for the Oo = 0.7 models, which 
we regard as the least interesting models presented in this paper. 



23 



2 I \ I I \ \ I I \ \ \ I I \ \ I I \ \ I I \ 1 

' .2 .4 .6 .8 1 

Figure 12: Magnification probability Pm vs. Qq. Open circles: flat models (including Einstein- 
de Sitter); crosses: open models. 

5.3. Classification of tlie Images 

In Figure 13, we plotted, for illustration purpose, some of the images we have encountered. Each panel 
shows the image of a single circular source of angular diameter 1". We have plotted the location of the light 
rays on the image plane, and enlarged the dots representing these rays until the images look continuous. 
Individual rays can be seen along the edges of the images. In the absence of lensing, each image would be 
circular and would contain 190 light rays. 

Figure 13a shows the most common case: a single image, magnified and sheared. In such case, there 
is no galaxy along the line of sight. The magnification is caused by the background matter located near 
the beam, while the shear is mostly caused by the background matter and galaxies distant from the beam. 
Figures 13b and 13c illustrate the two most common cases of double image. Figure 13b shows a case of 
strong lensing. In such cases, the magnification is always large, in the range /i = 3 — 7, the brightness ratio 
between the bright image and the faint one is small, of order of a few, and the brightest image forms an arc, 
while the faintest image is more compact and invariably shows a spike that points toward the other image. 
This is a consequence of the particular density profile we assume for galaxies (eq. [5]). Figure 13c is a case of 
weak lensing. In such cases, the magnification is small, of order = 1.2—1.3, and the brightness ratio is very 
large. The effect of lensing is to "carve" a second, faint image out of the brighter one. Figure 13d shows an 
Einstein ring, caused by a lensing galaxy nearly aligned with the source. For such rings, the magnification is 
always very high, sometimes larger than 10. Figure 13e shows a very different, and unusual, Einstein ring. 
The ring is pinched at two different locations, where the width is only one light ray, and the magnification 
is quite small, less than 2. This case is actually similar to the double image case shown in Figure 13c, and 
one can go from one case to the other continuously by shifting the location of the source. Figures 13f-13h 
show cases of triple image. The images can form either a circular pattern (as in Fig. 13f), in which case the 
brightness ratios between the images are small, or a linear pattern (as in Figs. 13g and 13h), in which case 
one image is usually much brighter than the other two, that image being located either in the middle, as 
in Figure 13g, or on the edge, as in Figure 13h. The latter one is similar to some triple images generated 
by Makino & Tomita (1995, their Figs. 6b and 6d). Figures 13i and 13j show cases of quadruple image. 
These cases suggest that several galaxies are involved in the lensing. Figure 13k shows an interesting case of 
a double image, one of thcni being a ring. This is actually a combination of the cases shown in Figures 13b 
and 13d. There are two galaxies near the line of sight, one responsible for forming the ring, and the other 
responsible for forming the second image. Figure 131 shows a very rare case of an Einstein ring with two 
holes. This requires an alignment between the source and 2 intervening galaxies, each of them producing a 
hole. Figures 13m and 13n show even rarer cases of an Einstein ring with 2 holes and a secondary image, 
and an Einstein ring with 3 holes, respectively. We have encountered only one of each of these cases in all 
of our experiments. 
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Figure 13; Images of Icnsod circular sources, for a few interesting cases: (a) single image; (b) and 
(c) double image; (d) and (c) Einstein ring; (f), (g), and (h) triple image; (i) and (j) quadruple 
image; (k) Einstein ring with secondary image; (1) double Einstein ring; (m) double Einstein ring 
with secondary image; (n) triple Einstein ring; (o) Einstein ring with central image; (p), (q), and 
(r) complex images caused by multiple lensing. 



Figure 13o shows an Einstein ring with a central spot. The presence of this spot is a consequence of 
the particular model we use for the galaxies. We do not represent galaxies as singular isothermal spheres, 
but rather as nonsingular isothermal spheres with a central core. Such objects can produce three images if 
the angular separation between the object and the source is sufficient small (see SEE, pp. 244 and 396, and 
§5.5.1 below). In this case, the three images are located on the line going through the source and the lens 
on the celestial sphere, one image on each side of the lens, and the third image in the core of the lens. If the 
lens and the source are aligned, there are no preferred directions, and the two outer images will turn into a 
ring, while the third image, located in the core of the galaxy, will form a "spot." This explains the existence 
of images like the one displayed in Figure 13n, but it does not explain the existence of rings without spot, 
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such as the one displayed in Figure 13d. In our experiments, less than 1% of the rings have a eentral spot. 
This is simply a resolution effect. As SEF point out, for lenses modeled as nonsingular isothermal sphere, 
the third image, located in the core, is very faint. If the brightness of that image is less than one light ray, 
we cannot resolve it. Indeed, we often found rings containing a spot composed of less than five rays, which 
we regarded as being underresolved. 

Finally, in Figures 13p-13r, we plotted some of the strangest images we encountered. These complex 
images result from lensing by several galaxies, combined with shear caused by the background matter. 

The reader should keep in mind that most of these cases are extremely rare. 99.7% of the images fall in 
the category illustrated by Figure 13a, a single image, magnified and sheared. Of the remaining cases, the 
vast majority of them fall in the categories illustrated by Figures 13b 13d. The cases shown in Figures 13e- 
13r are all very rare. In §5 we will discuss in more detail the probability of occurrences of these various 
cases. 



5.4. The Shear Distributions 

There are several ways to compiitc the shear caused by gravitational lensing. The most direct way is to 
compute iteratively the magnification matrix along each light ray (SEF, §9.1.2), and compute its eigenvalues. 
This approach is not very practical for our experiments. We consider beams composed of light rays separated 
by 0.064". The magnification matrix would tell us the prop(Ttic« (rotation and deformation) of a source 
having that size. We are, however, considering sources of angular diameter 1", containing, in the absence of 
magnification, 190 light rays. We could always try to somehow combine the magnification matricies for all 
the rays inside an image to infer its properties, but there is a much simpler approach. Since we are actually 
resolving the shapes of individual images, we can estimate the shear simply by computing the aspect ratio 
of the images. 

To compute the aspect ratio of an image, we first compute the geometrical center of the image on the 
image plane, 

Vgc,i = JtYI ' (28) 
' k 



where r^ is the location of light ray k in image ?', and Ni is, as usual, the number of rays in image i. We 
then compute the 2-dimensional quadrupole tensor of the image, 



Q — ^ ^ (^fc.t ''gCjj) (''fc,i ''gc,j) 



A B 
B C 



(29) 



where the last equality defines the components A, B, C. The aspect ratio of the image is obtained from the 
eigenvalues Ai, A2 of Q, as follows, 

02 1[^ + C- [(^-C)2+4S2]V2/ • ^ > 

where oi and 02 are the long and short "axes" of the image. One interesting property of this expression 

is that in the case of an elliptical image, ai and a2 are the semi-major and semi-minor axes of the ellipse, 
respectively. The resulting shear distributions are plotted in Figures 14-17. Error bars have the same 
meaning as in Figures 6-9. 



5.4.1. The as Dependence 

Figure 14 shows the shear distributions for models with the same values of Qq, Aq, and Hq, and different 
values of ag. In the absence of lensing, the distributions would be ^- functions located at 01/02 = 1. As 
for the magnification distributions, the larger the departure from a J-function is, the stronger the effect 
of lensing is. The top panel shows various Einstein-de Sitter models with different values of as- As as 
increases, the peak of the distribution decreases while the high-tail of the distribution increases. This was 
expected, since the large-scale structure, whose amplitude is measured by as, is the primary origin of the 
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shear. Interestingly, the distributions are much less sensitive to the valu(^ of cts for the other models. For 
the fio = 0.2, Ao = model, the curves in the high-tail regions are very similar and cross each other several 
times, while for the flo = 0.2, Aq = 0.8 model the cases ag = 0.8 and cts = 1.0 are almost undistinguishable. 
Therefore, the dependence of the shear distribution upon (Ts is more pronounced for models with large CIq. 
We found a similar trend in Figure 6 for the magnification distributions. 




T \ r 



Qo = 0.2 (Tj = 0.6 2 




Figure 14: Distributions of aspect ratios for various combinations of fio, Aq, and H^^ showing 
the effect of varying erg. Error bars indicate the 1 — a uncertainty AP on the mean, for a few 
representative bins. See text for details. 

5.4.2. The Hq Dependence 

Figure 15 shows the shear distributions for models with the same values of Qq, Aq, and as, and different 
values of Hq. The curves in each panel are very similar, except for the case Hq = 55kms~^Mpc~^ in the 
top panel, for which the statistics are poor. We do not find any particular trend, and the ordering of the 
curves with Hq is not even monotonic. As for the magnification distributions shown in Figure 7, the absence 
of dependence upon Hq results from competing effects. With larger Hq, the mean background density is 
higher, increasing the effects of lensing, but the cosmological distances are shorter, decreasing the effects of 
lensing. 
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Figure 15: Distributions of aspect ratios for various combinations of fio, Aq, and ag, showing 
the effect of varying Hq. Error bars indicate the 1 — a uncertainty AP on the mean, for a few 
representative bins. See text for details. 

5.4.3. The SIq Dependence 

The top panel of Figure 16 shows the shear distributions for two matter-dominated models (Aq = 0) with 
the same values of Ho and ag, and different values of Clo- This situation is similar to the one encountered 
in §5.1.3 for the magnification distributions. All elements of gravitational lensing — cosmological distances, 
mean background density, large-scale structure — depend upon ilg, but the dependence upon the mean 
background density dominates, and consequently the distribution is wider (that is, lensing is stronger) for 
models with larger fl^. The bottom panel of Figure 16 shows comparisons between flat models (f2o + Ao = 1)- 
This is a totally different situation. As decreases, Aq increases, and the dependence of the cosmological 
distances upon Aq becomes the dominant effect, resulting in stronger lensing for models with smaller SIq. 

5.4.4. The Aq Dependence 

Figure 17 shows the shear distributions for two models with the same values of fio, Hq, and as, and 
different values of Aq. As Ao increases, the distributions become wider, indicating that the effect of lensing 
is stronger. A larger value of Ao results in larger cosmological distances, which is clearly the dominant effect 
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here. The; ciffcct of introducing a finite vahie Ag is very large for the Og = 0.2 model, but becomes weaker 
for larger values of fio, as the flat and open models become more similar. 
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Figure 16: Distributions of aspect ratios for various combinations of Aq, Hq, and ag, showing 
the effect of varying ^Iq. Error bars indicate the 1 — a uncertainty AP on the mean, for a few 
representative bins. See text for details. 

5.4.5. Large-Scale Structure versus Galaxies 

All results presented in this section are consistent with the ones presented in §5.1. The magnification 

distributions and shear distributions have the same, or similar, dependences upon the cosmological param- 
eters. This strongly suggests that magnification and shear are two different manifestations of one single 
physical phenomenon: weak lensing. In §5.1.5, we argued that the effect of galaxies on the magnification 
distributions is small, because galaxies cover only a small fraction of the lens planes. For the shear distri- 
butions, the situation is different. While magnification is primarily caused by the matter located near the 
beam, shear is primarily caused by the matter located away from the beam. In this case, the cross section 
of galaxies becomes irrelevant. The lack of contribution to the shear from galaxies comes from two different 
effects. First, even for the lowest mean background densities we consider, ^Iq = 0.2, the galaxies account for 
only 24% of the total mass. Second, as explained in great details in Paper I, each galaxy in the simulation 
resides on top of an extended, smooth, compensating "hole" of negative density, which is introduced to 
account for the matter that has been removed from the intergalactic medium to form this galaxy in the first 
place. This is necessary, otherwise the process of "adding" galaxies to the simulation would not conserve 
mass. As a result, the total mass galaxy -I- "hole" is zero, and galaxies have no influence beyond the 
radius of the hole, which we fixed at rhoie = 1 Mpc. Hence, only the (usually) small fraction of galaxies 
which are located within 1 Mpc of the beam can contribute to the shear. 



This is an improvement over the approach used by Jaroszyhski (1991), which involved an overall, uniform 
reduction of the mean background density. 



29 





1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 




a,/a2 

Figure 17: Distributions of aspect ratios for various combinations of Qq, Hq, and as, showing 
the effect of varying Aq. Error bars indicate the 1 — a uncertainty AP on the mean, for a few 
representative bins. See text for details. 

5.5. The Multiplicity of Images 

5.5.1. Practical Considerations 

For all cases of multiple images found in our experiments, we have computed the image separations, 
defined as the angular separation between the center of the images. In doing so, we had to decide which cases 
should actually be identified as multiple images. We have designed an algorithm with pattern recognition 
capability to automatically identify multiple images and compute their separations, as well as the brightness 
of each image. It turns out that in a very large number of cases, one of the images contains only one or two 
rays. In such cases, that image is demagnified by a factor of 100, and the brightness ratio between the two 
images is of order of 100 or more. It is extremely doubtful that an observer could resolve such a faint image. 
Indeed, as we will argue in §5.5.2 below, the vast majority of observed gravitational lenses probably contain 
an unresolved image. 

To keep the number of cases with multiple images at a manageable level, we excluded from this study 
all images containing less than 5 light rays. Therefore, a double image with one image containing less than 
5 rays is treated as a single image, a triple image is treated as a double imago, and so on. By imposing this 
restriction, we are postulating that an image with less than 5 rays could never bo resolved. However, we 
are not assuming that an image with 5 rays or more is resolvable. The 5-ray limit was introduced mostly 
for convenience, to facilitate the analysis, but as we will show in §5.5.2 below, the number of double images 
we predict is much too large to agree with observations, indicating that a realistic cutoff would have to be 
much larger than 5 rays. 

5.5.2. Theory, Simulations, and Observations 

We define P„ as the probability that a randomly selected source will have n images. In our experiments, 
each cell on the source plane is a potential location for a source, and cells are given equal probability. To 
compute Pn we simply divide the number of cases with n images by the total number of cells included in 
the analysis. We have performed a total of 3,798 experiments. For most experiments, 841 cells are included 
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in the analysis. In some cases, however, cells located near the edge of the beam must be rejected (see §3.3). 
Combining all experiments for all models, we have a total of 3,137,675 cells, each of them representing a 
potential source. Excluding any image that is composed of less than 5 rays, we found 10,728 double images, 
126 triple images, and 6 quadruple images. No cases with five images or more were found. This contradicts 
the odd-number theorem (SEF, p. 172, Theorem 1), as well as observations, which find N2 ^ N4 ^ N3 
(Kochanek et al. 1998). This is a consequence of limited resolution. In most realistic lens models, whenever 
multiple images are produced, one image is always very faint (see, e.g., SEF, pp. 58 and 175), often too faint 
to be observed. Furthermore, for lens models such as nonsingular isothermal spheres (the model used in this 
paper for galaxies), three images at most can be produced, and the faintest one is located in the core of the 
Icnsing galaxy, making it very difficult to observe. Our experiments do not have such limitations, but are 
limited by the number of light rays used. An unlensed source contains 190 light rays, and therefore images 
with a magnification fj, < 1/190 = 0.0053 are unresolved. 



5.5.3. Double Images 

For each model, we computed the probability P2 of finding a double image. The results are plotted 
in Figure 18. There arc 43 points in each panel, corresponding to the 43 different cosmological models 
considered. The values of P2 are quite high for the reasons explained above. There is a strong trend of P2 
to increase with Aq, as the top right panel of Figure 18 shows. 
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Figure 18: Double- image probability P2, versus ilo (top left panel), Aq (top right model), Hq 
(bottom left panel), and as (bottom right panel). Each point corresponds to one cosmological 
model. Symbols indicate the value of the cosmological constant; crosses: Ao = 0; open squares: 
Ao = 0.3; open circles: Ao = 0.5; filled circles: Ao = 0.8. 
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In order to study the variations of P2 with the other parameters at fixed Aq, we use different symbols to 
designate the different values of Aq. The large scatter in the values of P2 is clearly caused by the dependence 
of P2 upon Ao, as the various symbols indicate. On all panels, Aq = 0.8 models (solid circles) are concentrated 
at the top, while Aq = models (crosses) are concentrated at the bottom. Looking at the crosses in the top 
left panel, and looking separately at the various symbols in the bottom panels, we find no obvious trend, 
indicating that P2 is essentially independent of flo, Hq, and ag at fixed Aq. These results implies that (i) 
double images, and multiple images in general, are caused by galaxies and not by the background large-scale 
structure, (otherwise there would be a strong dependence upon as) and (ii) the strong dependence of P2 
upon Ao indicate that the cosmological distances are the dominant effect in multiple imaging. 

The most important and interesting property of double images, after their likelihood of occurrence, is 
the distribution of angular separations between the images. This will be the topic of §5.6 below. 

5.5.4. Triple Images 

We found 126 cases with triple images, out of a total of 3,137,675 potential sources. With such a small 
number of cases, we cannot make a quantitative determination of the probability P3 for occurrence of triple 
images, and its relationship with the cosmological parameters. Still, we have enough cases to study the 
general properties of triple images, and to identify trends. It turns out that several of the results found for 
the double images apply to triple images as well. 

We plotted the properties of triple images in Figure 19. The top left panel shows a histogram of the 
magnification distribution. The magnifications are quite high, usually in the range n = 2 — 7, with some 
cases having ^ > 10. Triple images are strong lensing events caused by massive galaxies, though the tidal 
field of the background matter and nearby galaxies might play a key role in determining the shape of the 
images. 

As for the double images, there is a strong trend for triple images to occur in models with a large 
cosmological constant. We studied 17 models with Ao ^ 0, and found cases of triple images in all of them 
but one: the model flo = 0.7, Ao = 0.3, Hq = 65kms~^Mpc~^, as = 0.9, for which we have performed 64 
experiments. We also studied 26 models with Ao = 0, and found cases of triple images in only 13 of them, 
in spite of the fact that we have performed more experiments for these models than the Xq models. The 
top right panel of Figure 19 shows the frequency of occurrence of triple images, defined as the number of 
cases (indicated by numbers) divided by the number of experiments, versus Aq. About one fifth of the cases 
(27 out of 126) were found in models with Aq = 0, but 71% of the experiments (2,699 out of 3,798) were 
performed with these models. The trend is clear: the frequency of triple images increases sharply with Aq. 
As for the double images, this indicate that the cosmological distances are the dominant effect. 

The bottom left panel of Figure 19 is the scatter plot of the brightness ratios B2/B3 and B1/B2, where 
Bi, B2, and S3 are the brightnesses of the brightest, intermediate, and faintest images, respectively. The 
dashed line corresponds to B2/B3 = B1/B2. The brightness ratios can be very large, up to ^ 100, and 
for this reason we decided to use logarithmic scales. Triple images located above the dashed line have 
B2/B3 > B1/B2, and are composed of "two bright images and a faint one," while triple images located 
below the dashed line have B2/B3 < B1/B2, and are composed of "one bright image and two faint ones." 
These latter cases outnumber the former ones 2 to 1. This shows that individual galaxies are not solely 
responsible for the formation of triple images, since a galaxy modeled as a nonsingular isothermal sphere 
could only produce a triple image composed of two bright images and one faint, possibly unresolved image. 
Tidal perturbation by nearby galaxies and by the background matter, and lensing by several galaxies, must 
play an important role. 

The bottom right panel of Figure 19 is a scatter plot of the minimum angular separation Smin and 
maximum angular separation Smax between the images, in arc seconds. We can identify three specific 
limiting cases: (i) In the circular limit, defined by Smin = Smax (top dashed line), the three images form a 
circular pattern, as in Figure 13f. (ii) In the linear limit, defined by Smin = Smax/2 (bottom dashed line), 
the three images form a linear pattern, with the two outside images equidistant from the central one. These 
arc; the kind of triple images that a nonsingular isothermal sphere would produce, (iii) In the hierarchical 
limit, defined by Smin Smax (horizontal axis), two images are very close to each other, with the third image 
located much farther. These would invariably correspond to cases of one bright image and two faint ones, 
since the two images located near each other would have to be faint in order not to overlap. As Figure 19 
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shows, none; of these hmits seems to be preferred. The; seatter is large, and the; only noticeable thing is 
that the strong hierarchical limit does not seem to occur very often, as very few points are located near the 
horizontal axis. 
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Figure 19: Properties of the triple images. Tof left panel: Histrogram of the magnification /i. Top 
right panel: Frequency of occurrence of triple images versus Aq. The numbers indicate the actual 
number of cases found. The frequency is that number divided by the number of experiments. 
Bottom left panel: brightness ratio (intermediate to faintest) versus brightness ratio B1/B2 

(brightest to intermediate). The dashed line corresponds to B2/B3 = Bi/B2- Bottom right panel: 
minimum angular separation Smin versus maximum angular separation Smax- The dashed lines 
correspond to the "circular limit" Smin = Smax and the "linear limit" Smin = S'inax/2. 



5.5.5. Quadruple Images 

Quadruple images are extremely rare. We only found 6 cases, and two of them are plotted in Figure 13. 
All cases were found in models with Aq = 0.8, reinforcing the trend of multiple images being more frequent 
in models with large cosmological constant. The magnifications are jj, = 2.119, 2.910, 4.846, 8.344, 11.222, 
and 14.012. With 3 cases out of 6 having fx > 8, quadruple images are even stronger lensing events than 
triple images. According to the odd-number theorem, each of these cases must be at least a quintuple image 
with one faint, unresolved image. This is most likely the case for the quadruple image shown in Figure 13j. 
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These images forms a circular pattern, and there is probably a fifth, unresolved image in the center. More 
complex cases, such as the one shown in Figure 13i, could very well have more than 5 images. 

5.6. The Distribution of Image Separations 

5.6.1. Histograms of the Image Separations 

For each case with multiple images found in the experiments, we computed the angular separation 
between the images. The center of each image is computed using equation (28), and separations are computed 
between image centers. Figures 20 and 21 show histograms of the angular separations in arc seconds, for 
all models. All histograms are plotted on the same scale to facilitate comparison. On each panel, we have 
indicated the total number N2 of double images found. This number varies among the different models, 
first because some models are more likely to produce multiple images than others, and second because the 
number of experiments performed was not the same for all models. Models with A^2 < 200 tend to have 
very noisy distributions, and many more experiments would be required in order to determine the precise 
shape of these distributions. Still, several trends are apparent. We are considering sources with an angular 
diameter of 1". Since one image at least must be brighter than the unlensed source would be (SEF, p. 172, 
Theorem 2), the smallest possible image separation is 0.5" corrc^sponding to a double image like the one 
shown in Figure 13c. Most histograms in Figures 20 and 21 show a distributions that rises sharply from 
0.5" to 1", and then drops slowly at larger separations, with a high-tail that extends to separations of order 
4" - 6". 

As in the case of the double-image probability P2 discussed in §5.5.3 above, we find no obvious correlation 
between the shape of the histograms and the value of ug. This again indicates that double images are caused 
primarily by direct interaction between the beam and individual galaxies, and not by the large-scale structure. 
There is, however, a relationship between the largest angular separations and the value of Aq. For models 
with Aq = 0, the high-tail of the distribution function rarely extends beyond 4", while for Aq = 0.8 models, 
the high-tail often extends to separations of 6". As for the probability P2, the shape of the high-tail depends 
strongly upon the cosmological distances. Increasing these distances results in higher image separations 
for a given lensing galaxy. This affects the magnification distribution, by extending the high-tail to higher 
separation, and also the probability P2, by "separating" images that otherwise would have overlapped and 
been detected as a single image. Fukugita et al. (1990) found that the image separations were independent 
of the cosmological constant, but only considered the rms value {s^y^^ of the separation, and not the actual 
distributions. Since cases with large separations are rare, the details of the high-tail can hardly affect the 
rms value, and therefore these results are not in contradiction. 

In several histograms, especially the ones for Aq > models, we see a secondary peak at large separation. 
Consider for instance the model flo = 0.2, Aq = 0.8, Hq = 65 kms~''^Mpc~''^, erg = 0.8, which is indicated by 
an asterisk in Figure 20. There are no double images with separations between 4.00" and 4.75", but there 
are several images with separations larger that 4.75". This might seem like a very small effect that could 
be dismissed as a statistical fluctuation, but this feature is found in many histograms, suggesting that it 
could actually be real. This could possibly result from a coupling between galaxies and large-scale structure. 
Galaxies are predominantly responsible for multiple imaging. But most galax;ies are located inside clusters, 
where the density of background matter is high. This background matter might amplify the lensing effect 
of the galaxy, resulting in a peak at high separation angles. This issue requires more investigation. 

5.6.2 A Synthetic Angular Separation Distribution 

To gain more insight into the origin and properties of the angular separation distributions shown in 
Figures 20 and 21, we will now synthesize an angular separation distribution using a simple analytical 
model. We consider the fio = 0.2, Aq = 0.8, Hq = 65kms^^Mpc~^, erg = 0.8 model, and make the following 
assumptions: (1) Lensing is entirely caused by galaxies; we ignore the lensing by the background matter. 
(2) Each galaxy acts as if it was alone; we ignore the tidal effects of nearby galaxies, and the possibility 
of lensing events involving several galaxies. With these assumptions, the problem is reduced to studying 
lensing by isolated, nonsingular isothermal spheres. A problem of such importance has generated a great 
deal of interest in the past, and several analytical results have been derived (Dyer 1984; Hinshaw & Krauss 
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1987; Blandford & Kochanek 1987; Kochanek & Blandford 1987) which we can now apply to our model. We 
use the notation of SEF, §12.2.3. 




Figure 20: Histograms of the distribution of image separations in arc seconds. The value of the 
cosmological parameters and the number of double images are indicated in each panel (with Hq 
in units of kms~^Mpc~^). The large asterisk indicates the distribution which is replotted in the 
bottom panel of Fig. 22. 
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Figure 21: Same as Fig. 20, for different cosmological models. 

For each galaxy, we introduce a length scale ^o, defined by 

(SEF, eq. [8.34a], with v = V^a^), where v is the circular velocity given by equation (8). We now define 
the scaled quantities Xc = t'c/^o and y = r]/^o, where rc is the core radius given by equation (6), and r] is 
the distance between the source and the optical axis going through the observer and the center of the lens. 
We also define, for < 1, a critical radius = (1 — Xc^^)^^'^- The nonsingular isothermal sphere has the 
following properties (SEF, §12.2.3): (1) If Xc > 1, the source will have only one image. (2) If Xc < 1 the 
source will have one image ii y > ijr, and 3 images if y < y^- Hence, each lens which satisfies the condition 
Xc < 1 has a cross section for multiple imaging equal to ny^. 

In our analytical model, we consider as potential lenses all the galaxies located between the source and 
the observer. 12% of the galaxies, located mostly at redshifts < 0.1 or 2; > 2, are rejected as unable to 
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produce niTiltipk^ images, having scaled core; radii Xc > 1. The top panel of Figure 22 shows the distribution 
Xc for the remaining galaxies. The distribution is bimodal, with most early type galaxies (ellipticals and 
SO's) having Xc < 0.05, and most spiral galaxies having 0.2 < Xc < 1. In our model, spiral galaxies have core 
radii that are typically 10 times larger and circular velocities that are typically twice smaller than early type 
galaxies, as Table 1 shows. Since Xc oc rc/v^, we expect the values of Xc to be typically 40 times smaller 
for early type galaxies than for spiral galax;ies, explaining the bimodal distribution shown in Figure 22. The 
spread in the distribution within each peak is caused by the dependence of Vc and v upon the luminosity L. 
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Figure 22: Synthetic angular separation distribution for flo = 0.2, Ao ~ 0.8, Ha = 65kms~^Mpc~^, 
as = 0.8 model. Top panel: distribution of scaled core radii Xc- Middle panel: distribution of 
critical radii yr- Bottom panel: synthetic angular separation distribution (solid line), and actual 
distribution, taken from the thirteenth panel of Fig. 20 (dotted line). See text for details. 

The second panel of Figure 22 shows the distribution of critical radii. Small values of Xc lead to large 
values of yr, and vice versa. Early type galaxies have critical radii concentrated near the maximum possible 
value yr = 1,^^ while spiral galaxies have yr < 0.55. Hence the cross section for multiple imaging is dominated 
by early type galaxies. 

If a galaxy modeled as a nonsingular isothermal sphere produces multiple images (that is, 3 images), 
the angular separation s l)etw(H;n the two outermost images depends upon the source position y, and has a 
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//. 1 would require rc = 0, that is, a singular isothermal sphere 
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maximum value given by 

s^^ = ^{l-xiy/\ (32) 

Hinshaw & Krauss (1987), and Cheng & Krauss (1999) showed that the dependence of s on y is weak. 
Following the suggestion made by SEF (p. 396), we will assume that whenever multiple images occur, the 
image separation is of order s ~ Smax- We can then compute the distribution of image separations directly 
from the distribution of scaled core radii Xc- In doing so, we found out that many separations are significantly 
smaller than the source size, in which case the images strongly overlap and would be observed as one single 
image. To build a realistic distribution, we must impose limits on the smallest possible image separation 
that allows individual images to be resolved. We assume that sources have an angular diameter of 1". The 
smallest possible image separation is therefore 0.5", corresponding to an image configuration such as the one 
shown in Figure 13c. However, with such small image separations, the images will often overlap (Fig. 13c 
is a particular, "lucky" case). We assume that at separations s < 0.5" the individual images can never be 
resolved, that at separations s > 1" they can always be resolved, and that at separations 0.5" < s < 1" 
they can sometimes be resolved, with a probability that varies linearly from to 1 between s = 0.5" and 
s = 1". The requirement that s must exceed 0.5" in order to possibly create resolvable multiple images 
eliminates as potential lenses about half of the early time galaxies, and 97% of the spiral galaxies. Hence, 
early types galaxies are much more likely to produce resolved multiple images than spiral galaxies. Not only 
the early type galaxies capable of producing resolvable multiple images outnumber the spiral galaxies 20 to 
1, but their cross sections for multiple imaging are larger, as the second panel of Figure 22 shows. Notice, 
however, that both the numerical simulations and the analytical model treat galaxies as isothermal spheres. 
Blain, Moller, & Mailer (1999) have shown that the presence of a galactic disc can significantly increase the 
cross section for multiple imaging. Hence, our simulations and analytical model probably underestimate the 
contribution of spiral galaxies to multiple imaging. 

We computed the distribution of angular separations, using the above criterion for "resolvability." We 
also gave to each galaxy a weight w = to take the effect of cross section into account. The resulting 
distribution is shown in the bottom panel of Figure 22, by the solid line. For comparison, we plotted 
the actual distribution (dotted line) for this particular model, taken directly from the thirteenth panel 
of Figure 20 (indicated by an asterisk). The agreement is quite remarkable. The synthetic distribution 
reproduces the main features of the actual distribution: a sharp rise a separations s < 1", and a slow decline 
at separations s > 1". The analytical model underpredicts the number of cases with separations s ^ 1" — 1.6", 
and overpredicts the number of cases with separations s > 2.4". This is probably a consequence of ignoring 
the presence of the background matter. 

For the particular cosmological model considered in this section, wc found 787 cases with double image. 
Because this number is quite high, we are confident that the distribution of image separations is fairly well 
determined by the experiments. Hence, the bottom panel of Figure 22 serves to validate the analytical 
model. Armed with the knowledge that the analytical model can reproduce, to a reasonable accuracy, the 
actual distribution of separations, we intend to performed a much more detailed study of these distributions 
and their relationship to the cosmological parameters, based on the analytical model instead of ray tracing 
experiments. This study will be presented in a forthcoming paper (Martel, Premadi, & Matzner 2001). 

5.7. Einstein Rings 

Einstein rings are very common in our experiments, but most ot them are rather unspectacular: The 
image contains a few hundred rays, and the hole is made of a small number of "missing" rays, often less than 
10. Such holes have an angular diameter Dhoie of order 0.1" or less, and it is doubtful that observations 
could possibly resolve such small holes. Only the most massive galaxies can produce spectacular rings with 
hole angular diameters of order 1", such as the one shown in Figure 13d, but since only a small fraction of 
galaxies are very massive, such spectacular rings are quite rare. 

We need to decide where to draw the line between holes that are too small to be resolved, and holes 
that are not. This is obviously a very subjective decision, especially since the possibility of resolving small 
holes depends upon the details of observation and instrumentation. We decided, quite arbitrarily, to exclude 
from this study rings with a hole diameter Dhole < 0.64". Because this choice is arbitrary, we cannot 
derive precise statistics of the occurrence of rings and its dependence upon the cosmological parameters. We 
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can only describe the intrinsic properties of the rings, and express their relationship with the cosmological 
parameters in terms of general trends. 

The properties of the Einstein rings are summarized in Figure 23. The top left panel shows a histogram of 
the number of rings found versus their magnification fi. This shows that rings are usually high-magnification 
events. The histogram peaks at a magnification 5, and then drops slowly with magnification. Several 
rings have a magnification > 10, and the largest magnification we encountered was = 24.097 (it is the 
double ring shown in Fig. 13m). This justifies our use of the term "spectacular" to designate such rings. 
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Figure 23: Properties of the rings. Top left panel: histogram of the number of rings versus magni- 
fication jji. Top right panel: schematic plot showing the image of an unlensed source and a "typical 
ring," obtained by averaging the properties of all the rings found in the simulations. The horizontal 
bar indicates 1 arc second. Bottom panel: scatter diagram of the hole diameter Dhoie versus ring 
diameter Dring mcasiu-ed in arc second. The different symbols indicate models with and without a 
cosmological constant. The solid curves, form top to bottom, correspond to magnifications fJ. = 2, 
4, 8, and 16. The dotted lines, from left to right, correspond to thickness parameters Tring = 0.00, 
0.25, 0.50, and 0.75. 
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The top right panel of Figure 23 shows a "typical" ring, obtained by averaging the ring diameters Dring 
(the outer diameter of the images) and the hole diameters -Dhole over all rings found in all experiments. The 
average hole diameter is (Dhoie) = 1-09", slightly larger than an unlensed source. The average ring diameter 
is (Dring) = 2.62". 

In the bottom panel of Figure 23, we plot the distribution of ring diameters and hole diameters in arc 
seconds. The solid curves are contours of constant magnification, from left to right: = 2, 4, 8, and 16. 
Most rings have magnifications between = 4 and /i = 8 as we saw also in the histogram. The circles and 
crosses indicate models with Aq = and Aq ^ 0, respectively. Rings with /U < 8 are found in comparable 
numbers in both cases, but rings with > 8 are found predominantly in models with Aq ^ 0, and particularly 
in models with Xq = 0.8. This is the most important trend we found. The largest magnifications we have 
encountered, /z = 19.446 and /z = 24.097, were for models with Aq = 0.8. We also found a weaker trend: /U 
tends to increase with Hq. The dashed lines are, from top to bottom, lines of constant thickness parameter 
T'ring = 0.00, 0.25, 0.50, and 0.75, where the thickness parameter is defined by 

rp _ Dring ~ Dhole 

-tring — ^ • 1,<3<3J 

-^ring 

An infinitely thin ring has i?hoie = -Dring and Tring = 0; an infinitely thick ring has -Dhole = and Tring = 1- 
Figure 23 shows that most rings are quite thick, having T^ing in the range 0.5 — 0.75. The typical ring shown 
in the top right panel has a thickness parameter T^ing = 0.58. 



6. DISCUSSION 

The algorithms used for generating the density structures in the universe and for computing the trajec- 
tories of light rays are based on a certain number of approximations. In this section, we review these various 
approximation, and their possible effect on the results. 

The #1 approximation of any multiple lens-plane algorithm is that a continuous distribution of matter 
can be approximated by a finite number of lens planes, and that the results converge if the number of planes 
is sufficiently high. This issue was investigated by Lee & Paczynski (1990), and more recently by Martel et 
al. (2000). This latter paper showed that for sources located at redshift 2: = 3, 14 lens planes are sufficient 
to achieve convergence. In this paper, we used between 26 and 69 planes, depending on the model. Hence, 
the convergence of our multiple lens-plane algorithm is not an issue. 

Our algorithm for locating galaxies inside dense regions, which is described in detail in Paper I, ignores 
the existence of small-scale correlations, which lead to galaxy pairs such as the Local Group. However, as 
Figure 10 (bottom panel) shows, pairs of nearby galaxies are actually quite common in our simulations. Most 
of these pairs are not real physical pairs, being made of galaxies at di£Fc!rc!nt distances seen in projection, 
but this does not matter since only the projected surface density enters in the lens equation. The projected 
distributions of galaxies would hardly look any different if pairs were explicitly included in the algorithm. 

The relation given by equation (8), relating the circular velocity to the luminosity, might break down 
in the limit of small luminosities. This would not affect the magnification and shear distributions, which are 
determined primarily by the background matter, but it could affect the distribution of image separations. 
However, using the equations given in §5.6.2, we can easily show that the high-luminosity galaxies have larger 
cross sections for multiple imaging, and produce larger separations. With a Schechter luminosity function, 
low-luminosity galaxies are exponentially more abundant, but this does not affect the results because these 
galaxies produce separations that are much too small to be resolved. As we indicate in §5.6.2, imposing a 
minimum separation of 0.5" "eliminates" half of the early type galaxies and 97% of the spiral ones, and most 
of these galaxies are low-luminosity. 

Galaxies evolve. Not only their number density changes as a result of merging, but their individual 
density profiles might change also. We decided to ignore these effects in the current paper. Our method for 
locating galaxies in the computational volume is only approximate, and one might question the relevance of 
adding detailed modeling of the merging history and galaxy evolution, considering the uncertainties already 
present in the algorithm. We believe that our neglect of galaxy evolution and merging does not affect our 
results significantly. First, the weak lensing (magnification and shear) is unaffected by the galaxies. As for 
the strong lensing, the effect can be important. However, one critical point is that this paper presents a 
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com,parative study of cosinological models. Some of the lensing properties might change if the algorithm is 
modified. But these changes would not necessarily affect the trends seen in the variations of these properties 
with the cosmological parameters. To give a specific example: consider the double image probability P2 
plotted in Figure 18. Would the values of P2 be different if we had used different models galaxies? Most 
certainly. Would the trends revealed by Figure 18 {P2 increasing strongly with Aq, with no dependence 
on flo, Ho and ag) be different? Very unlikely. The focus of this paper was to determine the relationship 
between the lensing properties and the cosmological parameters. We intend to present a study of the effect 
of galaxy merging, and of considering various radial density profiles of galaxies, in a forthcoming paper. The 
focus of that paper will be different, and we will certainly not consider 43 different cosmological models. 

We superpose on top of each galaxy a "hole" of negative density to represent the matter that has been 
removed from the intergalactic medium by the galaxy formation process (Paper I). This is certainly the most 
uncertain part of our algorithm. This uncertainty simply reflects our lack of understanding of the galaxy 
formation process. At large scales, the density structures are dominated by the distribution of background 
matter, and the formation of these structures in CDM universes is well understood and accurately simulated 
by the P^M algorithm. At small scale, the density structures are dominated by galactic halos, and we 
represent these halos using analytical models that are in good agreement with observations, at least at low 
redshifts. It is the intermediate scale, immediately above the galactic scale, which is poorly understood. 
Galaxies must clearly remove matter from the intergalactic medium when they form, but the effect of this 
removal on the density structure of the intergalactic medium aroimd galaxies is not known. Yet. it cannot 
be ignored. Adding galaxies to the system using a Monte-Carlo method increases the total mass of the 
system, and therefore the mass of the background matter must be reduced somehow in order to conserve 
mass. By siiperposing a hole on top of each galaxy, we assume that when a galaxy forms, the mass that 
accumulates in that galaxy comes from its vicinity. This is highly uncertain. Fortunately, the effect of the 
holes on the lensing properties of galaxies is small. When the beam hits a galaxy, it encounters 3 different 
density structures: the galaxy itself, the hole superposed on top of the galaxy, and the background matter. 
Of all these three components, the hole is by far the least important one. The (negative) mass of the hole 
is equal to the mass of the galaxy, which is much smaller than the mass of the nearby background matter, 
even for models with CIq — 0.2. Furthermore, the holes are "spread" over radii much larger than the galaxies 
themselves. These holes should be seen as a small correction to the density of the background, to account 
for the removal of matter by the galaxy formation process. Notice also that according to the recent paper 
by Cheng & Krauss (1999), the effect of the background matter on the image separation caused by galaxies 
is at most 10%, and the presence of the hole introduces only a small correction to this value. 

We assume throughout this paper that sources have an angular diameter of 1". To test the effect of 
the source size on the results, we selected a particular model, flo = 0.2, Aq = 0.8, Hq — 65 kms~^Mpc~^, 
as = 0.8, and performed an additional 101 experiments using a smaller beam, 10.95" x 10.95" instead of 
21.9" X 21.9". With the same number of cells in the beam, this amounts to considering sources of angular 
diameter 0.5". We found that the magnification distribution P(/u) is independent of the source size, the 
distributions for 1" sources and 0.5" sources being nearly identical. Hovewer, the image multiplicities are 
reduced. Prom (^2,^3,^4) = (0.00882,0.00009,0.00001) for 1" sources, they dropped to (P2,P3,f4) = 
(0.00616, 0.00006, 0.00000) for 0.5" sources. This was expected, since a smaller source does not sample as 
well the inhomogeneities in the density structure which are responsible for lensing. 

7. SUMMARY AND CONCLUSION 

We have studied the propagation of light in inhomogeneous universes, for 43 different COBE-normalized 
CDM models with various combinations of f2o> Aq, Hq, and erg. We have performed a total of 3,798 numerical 
experiments, using a multiple lens-plane algorithm that enabled us to study the properties of weak and 
strong lensing simultaneously. Each experiments consisted of propagating a square beam of angular size 
21.9" X 21.9", composed of 341 x 341 light rays, from the observer to a source plane located at redshift z = 3. 
Our main results are the following: 

(1) At fixed Qo, Ao, and Hq, the magnification distribution depends upon as- As as increases, the 
low-tail of the magnification distribution shifts toward lower magnifications, because light rays are more 
likely to propagate through mostly underdense regions. The high-tail of the magnification distribution is 
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hardly affected. This rcsuh indicates that it is the background matter, and not galaxies, that are primarily 
responsible for the magnification of sources. 

(2) At fixed flo, Ao, and as, the magnification distribution becomes narrower as Hq increases, reflecting 
the fact that cosmological distances become shorter. This trend is rather weak, because the distance effect 
is partly compensated by the fact that the mean background density increases with Hq. 

(3) At fixed Ho and as, the low-tail of the magnification distribution shifts to lower values (that is, 
larger demagnification) as fio increases for Ao = models, and shifts to higher values as fio increases for 
^^0 + Ao = 1 models. In the case Ao = 0, the dependence of the mean background density upon Qq dominates 
over the respective dependences of the cosmological distances and large-scale structure upon Qq. In the case 
f^o + Ao — 1, the dominant effect is the dependence of the cosmological distances upon Aq. 

(4) At fixed Qq, Hq, and as, the magnification distribution becomes wider (that is, stronger Icnsing) as 
Ao increases, because of the increase in the cosmological distances. The effect is particularly large for models 
with Ao = 0.8. 

(5) The magnification probability P„j is almost independent of as, for any combination of ^q, Aq, and 
Hq, indicating that does not depend strongly on the amount of large-scale structure. Our interpretation 
is that the beam travels through several underdense and overdense regions whose effects mostly cancel out. 
Increasing as makes the underdense regions more underdense and overdense regions more overdense, but 
their effects still mostly cancel out. The only trend we found is that Pm increases with Qq at fixed Ao, Hq, 
and as- This suggest that Pm depends primarily on the mean background density. 

(6) The shear distribution has essentially the same dependences upon the cosmological parameters as 
the magnification distribution. At fixed flo, Xq, and Hq, the shear distribution becomes wider with increasing 
as- At fixed flo, Aq, and as, the shear distributions is very unsensitive to the value of Hq- At fixed Hq 
and as, the shear distribution becomes wider with increasing Q,q for Ao = models, and and narrower with 
increasing fio (and correspondingly decreasing Ao) for f2o -|- Ao = 1 models. At fixed VIq, Hq, and as, the 
shear distribution becomes wider as Aq increases. 

(7) The similarities found between the properties of the magnification distribution and shear distribution 
suggests that both phenomena have the same origin: weak lensing. The dependence of the magnification 
distribution upon as is a clear indication that the large-scale structure in the background matter, and not 
individual galaxies, are responsible for determining this distribution. This occurs because galaxies cover only 
a small fraction of each lens plane, and are not too likely to be hit by the beam. This argument does not 
apply to the shear distribution, however, since shear is not caused by the matter being directly hit by the 
beam, but rather by the distant matter. In this case, the absence of significant contribution from galaxies 
originate from their small mass (at most 24% of the total mass is in galaxies), and by the presence, in the 
algorithm, of a compensating underdensity on top of each galaxy to account for the mass removed from the 
background matter during the galaxy formation process. The strong dependence of the magnification and 
shear distribution upon Ao indicate that while the large-scale structure is responsible for these effects, the 
magnitude of these effects depends strongly upon the cosmological distances. 

(8) The double-image probability P2 increases strongly with Aq. Wc found no clear dependence upon CIq, 
Hq, and as- The absence of dependence upon as indicates that individual galaxies, and not the background 
matter are responsible for forming double images, which constitute examples of strong lensing. The strong 
dependence upon Ao indicates that, again, the dominant effect is the cosmological distances. 

(9) The distribution of image separations has properties similar to the probability P2: a strong de- 
pendence upon Ao and no dependence upon as- The distribution rises sharply from 0.5", half the angular 
diameter of the source, to 1", and then drops slowly, down to separations of order 4" for Ao = models 
and 6" for Ao = 0.8 models. In many cases, we also found a small, secondary peak is the distribution at 
separations of order 5". 

(10) Using an analytical model, which assumes that multiple images are entirely caused by galaxies, we 
have generated a synthetic distribution of image separations for one of the model. This synthetic distribu- 
tion reproduces the actual distribution obtained from the simulations remarkably well, indicating that the 
assumption of galaxies being responsible for generating multiple images is correct. The only discrepancies 
are at separations of 1" — 1.6", where the synthetic distribution is too low, and at separation larger than 
2.4", where the synthetic distribution does not drop fast enough. These discrepancies probably result from 
the presence of the background matter, which is ignored in the analytical model. 
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(11) We only found 126 cases of triple images and 6 cases of quadruple images. With such small 
numbers, we cannot do any precise determination of the probabilities P3 and P4, but we can identify several 
trends. Triple and quadruple images are strong lensing, high-magnification events caused by galaxies. Their 
magnification distribution peaks around /i = 4, and the high-tail extends to values larger than fi = 10. 
These cases are predominantly found in Aq > models, indicating again that the effect of the cosmological 
distances dominates. Triple images come out in a variety of pattern, including circular and linear patterns. 
About 2/3 of the triple images are made of one bright image and two faint ones. 

(12) It is difficult to determine precise statistics for the occurrence of Einstein rings. Rings are very 
common, but we rejected most of them because the "hole" was very small. We only included in the analysis 
rings with hole diameters larger than 0.64". Rings are the most extreme cases of strong lensing found in 
our experiments. Most rings have magnifications in the range = 4 — 8, and the magnification distribution 
extends to /x = 25. The brightest image found in any given experiment is almost always a ring. Rings are 
usually quite thick, with a thickness comparable to the angular size of the source. Rings are found both in 
Ao = models and Aq > models, but bright rings, with /x > 8, are found predominantly in Aq > models. 

We can summarize these results as follows: (1) The cosmological distances play a critical role in nearly 
every aspect of gravitational lensing, both weak and strong. Consequently, the properties of gravitational 
lenses depend much more strongly upon the cosmological constant Ao than any other cosmological parameter. 
(2) Magnification and shear are examples of weak lensing caused primarily by the distribution of background 
matter, with negligible contribution form galaxies. Consequently, these effect are sensitive to the value of the 
rms density fluctuation as- (3) Multiple images and rings are examples of strong lensing, caused by direct 
interaction with galaxies, with at most a small contribution from the background matter. Consequently, the 
properties of multiple images are independent of as- They are determined by the cosmological distances, 
which depend primarily upon Aq, and by the details of the galactic models, which are usually independent 
of the cosmological parameters. Therefore, observations of weak lensing can be used to determine the 
cosmological constant and the unbiased density structure of the universe (that is, without having to assume 
some biasing factor between luminous and dark matter), while observations of strong lensing can be used to 
determine the cosmological constant and the internal structure of galaxies and clusters. (4) The dependences 
upon Hq and Slg are not as simple, because varying these parameters affects gravitational lensing in several 
ways that often partly cancel each other. For instance, the small dependence of the magnification and shear 
distribution upon Hq results from the competing effects of increasing cosmological distance while reducing 
the mean background density. Determining Ao and as from observations seems much more promising than 
determining Oo and Hq- 

Our experiments consider only sources located at redshift z = 3, in which case most of the matter 
responsible for lensing is located at redshifts z ^ I. This is not an important limitation of our study. As 
Figure 5 shows, the angular diameter distances vary only weakly with redshift in the limit z > 1. Hence, 
our conclusions remain valid for sources located at larger redshifts. 

In conclusion, our study shows that it is difficult to single out the effect of each particular cosmological 
parameter. This supports the idea of conducting this simultaneous cosmological parameter survey. 
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